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Préface 


Le rayonnement électromagnétique des dispositifs électroniques modernes est une 
préoccupation majeure des concepteurs et utilisateurs, soucieux de respecter les 
règles de la compatibilité électromagnétique. 

Les mesures directes, bien que parfois délicates, ne présentent pas de problèmes 
insurmontables pour les laboratoires compétents au sens de la directive CEM. La 
difficulté concerne en réalité, les études de nature prédictive, qu’il conviendrait de 
conduire par simulation avant même le routage des cartes. Les circuits analogiques 
utilisés en traitement de signal, seraient alors examinés principalement dans le 
contexte de l’immunité ou de la vulnérabilité, alors que les cartes numériques et les 
cartes de puissance à découpage H.F., le seraient principalement sur le plan de 
l'émission, respectivement en champ électrique et magnétique. 

Ces études sont difficiles et ne sont actuellement développées que par quelques 
laboratoires de recherche. Elles nécessitent des moyens informatiques puissants et 
un investissement de temps considérable. 

J.L. Cocquerelle et C. Pasquier, dans le cadre du présent ouvrage, présentent une 
version didactique, donc accessible sur le plan pédagogique, des résultats de travaux 
de recherche autorisant une approche simplifiée du rayonnement des cartes de 
puissance à découpage. Ils ont en effet créé dans ce but un concept original, celui 
des «flots rayonnants », qui correspondent sur une carte électronique à des zones 
émettrices de champ magnétique à contribution prépondérante. 

Un support informatique est joint au livre, il contient un grand nombre de 
représentations graphiques trés pertinentes permettant la meilleure compréhension 
possible au lecteur. 


Je tenais, en préfacant cet ouvrage, à souligner l'effort que ces deux chercheurs de 
l'université de Nantes ont entrepris pour assurer le transfert vers la pédagogie, de 
travaux de recherche dans un domaine édicté et soutenu par l'évolution des 
technologies en électronique. 


Joseph SAILLARD 


Vice-Président du Conseil Scientifique de l'Université de Nantes 
Chargé de la Recherche 


Note des auteurs 


L'évolution des technologies en électronique permet de réaliser des convertisseurs 
de puissance à fréquence de découpage de plus en plus élevée. L'étude des champs 
émis par de telles structures devient donc une obligation pour le concepteur soucieux 
d'éviter la pollution électromagnétique de l'environnement proche et lointain. 

Cet ouvrage résulte de travaux de recherche théorique et appliquée, que nous 
avons conduits dans cette perspective, au laboratoire d'électronique de puissance 
et de compatibilité électromagnétique (LEPCEM), de l'école polytechnique de 
l'Université de Nantes. 

La représentation des champs, dans un livre, simplement réalisé en monochrome, 
est forcément réductrice en qualité graphique. Nous avons donc conçu un support 
informatique, permettant de visualiser en couleurs et avec une meilleure définition, 
la majorité des illustrations du livre ainsi que de nombreuses autres venant en 
complément. 


Le mode d'emploi est trés simple : 


e L'icóne « E, », portée à droite d'une illustration du livre, indique que celle- 
ci est reprise sur le support informatique à la page x de celui-ci ; 


e L'icóne « ЫСЮх-у », portée à gauche dans le texte suivant une illustration du 
livre, indique que des représentations complémentaires sont disponibles sur le 
support informatique de la page x à la page y ; 


• L'icóne « LHI Page x », portée en bas et à droite d'une illustration du support 
informatique, indique que celle-ci est présente ou annoncée, dans le livre à la page x. 
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INTRODUCTION 


Cet ouvrage s’inscrit dans les activités de compatibilité électromagnétique (CEM) 
du génie électrique. Il concerne principalement les convertisseurs de puissance à 
fréquence de découpage élevée, que l’on trouve par exemple dans les alimentations 
des micro-ordinateurs, celles des matériels audio-vidéo modernes, et tous les 
dispositifs électroniques de manière générale. Ces structures sont intrinsèquement 
génératrices de perturbations conduites et rayonnées, ces dernières pouvant altérer le 
fonctionnement de circuits sensibles proches dont, parfois, la carte de commande de 
la structure elle-même (autoperturbation). 

Les études en CEM concernent généralement le rayonnement à une distance 
suffisante pour obtenir une onde plane et ne retenir alors que le champ électrique. 

Plusieurs d’entre elles s’intéressent plus particulièrement aux convertisseurs de 
puissance et à la prédiction de leur rayonnement. Les méthodes numériques 
employées (méthode des moments, des éléments finis...) ont souvent pour objectif 
la modélisation d'ensemble du convertisseur autorisant l'obtention principalement 
du champ électrique global émis. Ces études sont complexes et s'adressent à un 
public averti. 

Certains logiciels de CEM commercialisés, essentiellement destinés aux cartes 
numériques, permettent d'obtenir en 2D et parfois en 3D la cartographie du champ 
électrique résultant. Ces logiciels sont spécifiques et onéreux ; les moins coüteux 
d'entre eux procédent parfois à des simplifications extrémes basées sur l'estimation 
dans le « pire des cas » au sens des normes en CEM. 

Si l'on se préoccupe du rayonnement électromagnétique dans l'environnement 
immédiat d'une carte électronique numérique, on sait que le champ électrique sera 
prépondérant car les circuits utilisés fonctionnent avec des temps de commutation 
trés courts, (quelques ns) et les gradients de tension correspondants (dv/dt) peuvent 
excéder 10? V/s, alors que les courants étant trés faibles (unité courante : le WA) le 
champ magnétique proche est trés réduit. S'il s'agit en revanche, d'une carte 
électronique de puissance à découpage, les temps de commutation sont rarement 
inférieurs à 100 ns, mais, dés lors que les courants mis en jeu se chiffrent en 
ampères, le champ magnétique devient important, et son évaluation est nécessaire. 

Dans cette perspective, les auteurs de cet ouvrage ont mis à profit l'expérience 
pratique qu'ils ont acquise en laboratoire de tests en CEM concernant diverses cartes 
électroniques embarquées dans des convertisseurs de puissance à découpage. Cela 
leur a inspiré une démarche simplifiée particuliére. En effet, au voisinage immédiat 
de telles cartes, le champ magnétique mesurable avec une sonde adaptée est 
prépondérant, et se situe en des endroits particuliers qu'ils considérent comme des 
«îlots rayonnants ». 


Ces « flots » sont souvent repérables sur la sérigraphie du circuit imprimé 
correspondant. Les auteurs ont alors étudié la possibilité d’assimiler ces « îlots », au 
demeurant peu nombreux sur une carte standard, à de simples « boucles 
équivalentes ». Cette démarche apporte une simplification évidente accentuée par 
l'utilisation de calculs associés allégés ; elle peut sembler manquer de précision, 
mais apparaîtra plus pragmatique à la plupart des industriels concepteurs qui 
connaissent par expérience les « points chauds » de leur carte. 

Les auteurs ont donc conduit avec précaution une étude spécifique sur le 
rayonnement magnétique des boucles de puissance devant permettre d’obtenir un 
outil d'évaluation et de représentation du champ magnétique accessible à tous. 

Dès lors que le champ est connu en tout point de l’espace, il est possible de 
rechercher la force électromotrice qui pourrait étre induite dans une boucle sensible 
de l’environnement immédiat du convertisseur. En général, le concepteur ne 
souhaite pas mettre en place de blindage ni d'écran et dispose, le plus souvent, d'un 
coffret non métallique. Un apport qui a semblé original aux auteurs de ce livre dans 
ces conditions, concerne la détermination de l'orientation optimale de cette boucle 
sensible pour minimiser directement l'amplitude de la f.e.m. induite. 

Cela se traduit pour le concepteur de la facon suivante : 

- s’il se fixe un point de l’espace environnant, il disposera rapidement, après 

simple repérage et spécification des « îlots », des données de positionnement et 

d'orientation d’une carte sensible éventuelle. La perturbation sera alors minimisée. 

- s'il dispose de plusieurs emplacements disponibles à l’intérieur du coffret 
pour y fixer la carte sensible qu'il veut ajouter, il pourra obtenir une information 
sur les zones à retenir en priorité que nous appellerons par la suite « zones 
tranquilles ». Dans ces conditions, la f.e.m. induite dans une boucle victime sera 
déjà réduite, le choix d'une orientation optimale pouvant conduire jusqu'à 
pratiquement l'annulation. 


Le but de cet ouvrage est donc double : proposer au lecteur une approche simple 
concernant l'analyse des phénoménes électromagnétiques et les calculs associés, et 
le sensibiliser aux implications industrielles pouvant en découler. 


CHAPITRE I 


ASPECTS ESSENTIELS 
DE 
L'ÉLECTRONIQUE DE PUISSANCE 


Ce chapitre a pour but de présenter et de caractériser dans une démarche « CEM », 
les semi-conducteurs utilisés classiquement dans les convertisseurs de l'électronique 
de puissance. 

L'objectif final est de déterminer les formes d'ondes types des courants qui 
pourront étre à l'origine du rayonnement magnétique au sein d'un montage. 

Avant tout il est donc nécessaire de rappeler l'essentiel sur la commutation. 


I.1 La commutation en électronique de puissance 
Le fonctionnement des interrupteurs est déterminant au regard de l'importance des 
champs qui vont résulter de leur insertion dans un montage. 


1.1.1 Principes fondamentaux 

Les convertisseurs de l'électronique de puissance sont présents dans la quasi- 
totalité des appareils et systémes utilisant l'énergie électrique. 

Leur róle est de modifier la source électrique, c'est-à-dire le réseau distribué, soit 
alternatif, soit continu, pour satisfaire aux exigences édictées par la nature et les 
caractéristiques du récepteur. Les modifications peuvent concerner la nature méme 
de la source (courant ou tension) ainsi qu'un ou plusieurs de ses paramétres, dont 
essentiellement : l'amplitude, la fréquence, le régime de mise à la terre, le débit 
maximum. 

Le principe fondamental de l'électronique de puissance est de procéder à 
l'adaptation des sources uniquement par découpage temporel des tensions et 
courants, donc sans pertes appréciables d'énergie. 

Les composants, qui assurent les découpages, présentent alors deux phases de 
fonctionnement importantes au regard de la CEM : les commutations à l'ouverture et 
à la fermeture. Suivant les semi-conducteurs utilisés et la structure du convertisseur 
considéré, différents types de commutations sont à considérer : 

- concernant l'interrupteur, « la commutation » est le « mécanisme » physique 
d'ouverture ou de fermeture des semi-conducteurs constituant l’interrupteur. Elle 
est soit « commandée » c'est-à-dire provoquée par l'application d'un signal 
directement sur la commande du composant, soit « naturelle » ou « spontanée » si 
elle a lieu sans intervention extérieure. Dans ce dernier cas, elle résulte 
uniquement du passage par zéro de la tension appliquée au composant ou du 
courant qui le traverse ; 

- pour un convertisseur, le mot commutation est relatif aux événements 

concernant les éléments de la structure directement associés aux interrupteurs qui 

peuvent étre considérés non plus individuellement, mais assemblés dans des sous- 
ensembles appelés « cellules de commutation ». 

C'est au niveau de celles-ci que l'on introduit les notions de commutation « dure » 
et de commutations « douce » et « adoucie ». 
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1.1.2 Commutation dure 
La commutation dure, également appelée commutation forcée, résulte de l'emploi 


d'un composant commandé à l'ouverture et à la fermeture. Les contraintes 
importantes qui sont alors imposées à celui-ci dépendent fortement des inductances 
de connexion et de câblage et de diverses capacités parasites, ce qui entraîne : 

- des pertes élevées ; 

- des surtensions et des surintensités aux bornes de l'interrupteur ; 
de rapides variations de tension et courant (dv/dt et di/dt) ; 

- de fortes oscillations, en particulier pour le découpage à haute fréquence. 

Pour limiter ces effets, on peut associer aux semi-conducteurs des circuits d'aide à 
la commutation (CALC), dont le róle essentiel est de réduire les pertes 
correspondantes et de limiter les surtensions et/ou surintensités. 


1.1.3 Commutation douce 
En commutation naturelle ou douce, les interrupteurs présentent une commutation 
commandée et une commutation spontanée. 
Il est possible de regrouper ces convertisseurs en deux familles : 
= la première correspond à des commutations pilotées par les grandeurs d'état 
de nature électrique. Il s'agit par exemple de certains convertisseurs reliés au 
réseau alternatif. C'est également le cas des convertisseurs quasirésonnants qui 
utilisent un couple de composants passifs L-C pour provoquer l'annulation du 
courant ou de la tension. 
-  ]a seconde est relative aux structures pour lesquelles les commutations 
spontanées des interrupteurs sont obtenues consécutivement à la commutation 
d'autres interrupteurs, c'est le cas pour les hacheurs de puissance à un thyristor 
principal dont l'extinction est assurée par un circuit auxiliaire de nature oscillante 
comprenant un deuxiéme thyristor dit de soufflage. 


1.1.4 Commutation adoucie 

La commutation est adoucie par l'emploi de CALC associés aux composants, dont 
les deux mécanismes essentiels consistent à : 

- ralentir la montée en courant à la fermeture par adjonction en série d'une 

inductance avec l'interrupteur ; 

- ralentir la croissance de sa tension aux bornes par mise en parallèle d'un 

condensateur. 

Avec un CALC, le composant est à méme de commuter la totalité du courant alors 
que la tension à ses bornes reste encore faible ou, par dualité, supporter la totalité de 
la tension, le courant étant annulé. A noter pour les puristes qu'on ne peut 
théoriquement pas parler de commutation à zéro de courant (ZCS) ou zéro de 
tension (ZVS) car ces grandeurs ne font qu'évoluer différemment tout au long de la 
commutation et ne sont jamais totalement nulles. 


Chapitre I Aspects essentiels de l’électronique de puissance 


L'association de ces éléments passifs a bien entendu ses revers, comme le risque 
de surtension due à l’inductance à l’ouverture, et de surintensité due au condensateur 
à ja fermeture, ou encore celui de la génération d’oscillations indésirables. Un 
CALC est en conséquence un circuit plus complexe incluant d’autres éléments 
comme par exemple une résistance pour dissiper l'énergie liée aux oscillations ainsi 
que des diodes pour assurer les connexions unidirectionnelles indispensables. 


x 


La figure L1 présente un CALC dissipatif opérationnel à l'ouverture et la 


fermeture. 
Fri 


Figure I.1 : CALC dissipatif 


Pour les petites puissances, des CALC plus élaborés peuvent étre mis en ceuvre. 
Ils se composent généralement de condensateurs, d'inductances et d'un autre 
interrupteur. La commande de ce dernier permet de réaliser une oscillation locale 
qui provoque la commutation du composant principal au «zéro de tension » ou 
« zéro de courant ». Ce sont les techniques ZVS et ZCS évoquées précédemment. 

En conclusion, les CALC améliorent les commutations puisque les surtensions et 
surintensités sont réduites au niveau de l'interrupteur qui ne supporte plus seul les 
différents traumatismes. Cependant, on ne supprime pas les variations rapides de 
courant et de tension (di/dt et dv/dt) qui sont en fait reportées sur les composants 
réactifs des CALC, ce qui peut augmenter les probabilités de rayonnement 
indésirable, en particulier en champ magnétique, par les éléments bobinés. Il faut 
encore souligner que le renvoi éventuel à la source de l'énergie de commutation ainsi 
que celle des oscillations parasites peut réduire la qualité du signal en amont. 


I.2 Convertisseurs et composants associés 


Ce paragraphe est une présentation sommaire des différents convertisseurs, qui a 
pour objectif principal de rappeler quels sont les semi-conducteurs associés, et les 
perturbations essentielles qui les caractérisent. 


1.2.1 Convertisseurs à la fréquence réseau 

Il s'agit essentiellement des redresseurs commandés, des gradateurs et des cyclo- 
convertisseurs. 

Ces structures traditionnelles sont réalisées dans une trés large plage de puissance, 
elles utilisent essentiellement des thyristors commandés par la phase. Les 
commutations à l'ouverture des composants sont spontanées. Les vitesses de 
commutation sont peu élevées car limitées par la nature des sources et les 
caractéristiques des composants. 
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Les principales perturbations observées concernent essentiellement les 
harmoniques et sous-harmoniques des courants prélevés au réseau. Le rayonnement 
à proximité immédiate de ces structures peut être élevé, compte tenu des amplitudes 
très fortes des tensions et courants mis en jeu, mais la fréquence de découpage étant 
très basse, les champs décroissent très vite avec l'éloignement. On s'intéressera 
alors en CEM, essentiellement aux champs dits « proches ». 


1.2.2 Convertisseurs à fréquence industrielle (< 16 kHz) 


Il s’agit essentiellement des hacheurs de puissance utilisés par exemple dans les 
variateurs de vitesse pour moteurs à courant continu et des onduleurs destinés aux 
moteurs alternatifs, ou encore aux alimentations statiques de secours. 


1.2.2.1 Les hacheurs de puissance 


Les sources amont et aval sont de nature continue ce qui n'autorise pas bien sür la 
commutation spontanée des interrupteurs employés. Ceux-ci sont donc commandés 
à l'ouverture et à la fermeture. Les commutations générent alors des gradients de 
courant et de tension importants dans les différentes connexions du convertisseur. 

Les composants utilisés en trés forte puissance sont des thyristors avec circuit de 
soufflage, ou des GTO (thyristors blocables par la commande). En moyenne 
puissance, le transistor IGBT est préféré. 

Les perturbations conduites couvrent un spectre plus étendu que précédemment, 
puisque les fréquences de découpage choisies peuvent atteindre 16 kHz dans un 
souci de réduction des bruits acoustiques. 

Les vitesses de commutation sont encore peu élevées car toujours limitées par la 
nature des sources et les caractéristiques des composants. 

Les champs proches, du fait des amplitudes élevées des tensions et des courants, 
seront importants et avec un spectre un peu plus étendu en fréquence (pouvant 
atteindre quelques MHz) que pour les convertisseurs précédents. 


1.2.2.2 Les onduleurs 


Les onduleurs ont une structure en pont ou demi-pont et peuvent étre regroupés 
selon trois familles : les onduleurs de tension, les onduleurs de courant et les 
onduleurs a résonance (résonance série ou résonance paralléle). 

Les deux premières familles utilisent une commande à modulation de largeur 
d'impulsion (M.L.I.) unipolaire ou bipolaire ou génèrent des niveaux discrédités 
(marches d'escalier). Chaque commande a pour objectif la réduction des 
harmoniques de tension ou de courant aux basses fréquences. 

Les onduleurs à résonance sont sur charge réactive, ils utilisent un élément dual en 
série ou parallèle pour créer ainsi une source réceptrice oscillante, dont les passages 
à zéro des grandeurs électriques autoriseront les commutations naturelles. Les 
courants et tensions obtenus sont d’allures sinusoïdales. Les fréquences de 
fonctionnement doivent être proches de la fréquence d'oscillation du dipóle LC pour 
limiter leurs harmoniques. 
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Les onduleurs nécessitent généralement des interrupteurs bidirectionnels en 
tension ou en courant, entièrement commandés (GTO, IGBT ou transistors MOS 
associés à une diode). Les commutations peuvent avoir lieu plusieurs fois par 
période et être irréguliérement réparties (commande à M.L.L). C'est alors une 
séquence qui se répéte à la fréquence de fonctionnement. Il est ainsi difficile de 
parler de fréquence de commutation pour ces composants. On peut toutefois 
introduire la notion de « motif ». 

Les perturbations conduites sont observables sur l'entrée continue et sont 
directement liées aux découpage des ondes réalisant la sortie alternative. 

Le rayonnement à proximité dépendra de la puissance mais aussi du type des 
composants, puisque leur vitesse de commutation peut étre, comme nous le verrons 
ensuite, nettement supérieure dans le cas des MOS ou de certains IGBT de dernière 
génération. 


1.2.3 Convertisseurs à fréquence élevée (30 KHz - 100 kHz) 

L'utilisation de fréquences de fonctionnement élevées permet de s'affranchir des 
bruits acoustiques, de diminuer la taille des condensateurs, des inductances, des 
transformateurs et de faciliter le filtrage. Cependant, les composants capables de 
fonctionner à ces fréquences admettent des limites en courant et tension restreintes. 
Ces convertisseurs sont donc principalement utilisés en moyenne et petite puissance. 

Les alimentations à découpage constituent l'essentiel de cette famille, elles 
assurent principalement la conversion continu-continu et alternatif-continu (ajout 
d'un pont redresseur en téte). Techniquement, ces structures du type « hacheur », 
relient des sources de nature différente (hacheur sans accumulation) ou identique 
(hacheur à stockage intermédiaire). 

Les composants les plus employés sont les transistors MOS de puissance, les 
transistors bipolaires et les IGBT rapides. Les vitesses de commutation sont 
imposées soit par le composant soit par le CALC. En régle générale, les temps de 
commutation sont réduits au maximum pour conserver au convertisseur la meilleure 
dynamique possible et de ce fait le rayonnement électromagnétique peut étre 
important. Sa réduction est possible dans les cas particuliers des alimentations 
« quasirésonnantes » ou « résonnantes ». 


1.2.4 Caractéristiques des SC de puissance associés aux convertisseurs 

Dans cet ouvrage, nous nous limitons aux semi-conducteurs (SC) destinés aux 
petites et moyennes puissances et cherchons essentiellement les caractéristiques 
concernant leurs commutations et, plus particulièrement, celles qui déterminent les 
variations rapides du courant. 

Les caractéristiques fournies par les constructeurs sont variées, elles donnent 
généralement un aperçu des performances optimales du composant. Les tests sont en 
effet réalisés avec des circuits particuliers qui sont pour la plupart exempts des 
éléments parasites que l'on peut trouver sur des cartes industrielles (inductances de 
ligne, capacités entres pistes, éléments parasites d'autres composants...). 
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Compte tenu de cela, certaines valeurs relevables sur des montages industriels 
pourraient être beaucoup plus importantes que les niveaux proposés. 

Les caractéristiques présentées ci-après pour les différents composants signalés 
précédemment ont été prelevées dans différents recueils de constructeurs connus. 
Nous y avons sélectionné des semi-conducteurs classiquement utilisés et, pour 
chaque famille, des caractéristiques types, moyennées, en vue d’établir à terme une 
forme d'onde de référence du courant commuté. 


I.2.4.1 Le thyristor ou SCR 


Le thyristor (Silicon Controlled Rectifier) est un composant relativement lent, bien 
adapté pour des fonctionnements aux basses fréquences en forte puissance. C'est un 
composant unidirectionnel en courant et bidirectionnel en tension, dont la fermeture 
est commandée, alors que l'ouverture est spontanée au zéro de courant (figure I.2). 

K ir 


T GC 


Figure 1.2 : Thyristor 


Vr 


VT 


Le tableau L1 ci-dessous, accompagné de sa légende (figure 1.3), présente 
quelques caractéristiques sélectionnées de thyristors de différentes gammes. 


dv/dtd 


Tableau I.1 : Caractéristiques sélectionnées de quelques thyristors 
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Légende : 


Figure 1.3 : Ondes du thyristor 


-  t, (min) : temps minimum estimé pour atteindre le niveau 3 Irms(max) ; 

- а/а: vitesse maximale de croissance et de décroissance du courant ; 

- tą: temps de blocage du thyristor. C'est le temps minimum entre le passage à 
zéro du courant dans le sens direct et la montée de la tension directe sans remise 
en conduction du thyristor ; 

-  dv/dt: vitesse maximale de croissance de la tension directe du thyristor. 
Au-delà de cette valeur, il y a risque d’amorçage intempestif du thyristor. 


Nous proposons de retenir des ordres de grandeur qui permettront par la suite de 
définir l'onde type de courant. A noter que la puissance commutée n'atteint pas le 
maximum possible pour le thyristor, du fait que nous souhaitons pouvoir comparer 
les différents semi-conducteurs dans une plage de puissance commune. 

- temps de commutation à la fermeture : 0,1 us à quelques us ; 
-  di/dt: 30 A/us à 170 A/us (200 à 1 000 A/us en non répétitif). 


1.2.4.2 Le thyristor GTO 


Le GTO (Gate Turn Off), est trés utilisé dans le domaine des fortes puissances 
(figure I.4). La fermeture et l’ouverture sont commandées. 


K e ir 
i CG 


Figure 1.4 : GTO 


Quelques caractéristiques sélectionnées de thyristors GTO sont données dans le 
tableau I.2. 
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I (A) tu (US) | tmin (US) | di/dt ba (US) dv/dt 
max max 
(A/S) (Vius) 


CC 
ES 
CR 


Tableau 1.2 : Caractéristiques sélectionnées de quelques GTO 


Légende : 

= tg: temps de descente rapide du courant durant la conduction en inverse de la 
jonction gáchette — cathode ; 

= tą: temps nécessaire au courant de queue pour s'annuler (recombinaison des 
porteurs libres) ; 

— trmin# lLmwss€ 


di/dt max ' 


Le temps de commutation à l’ouverture peut valoir plusieurs microsecondes 
(try + (4), mais dans le contexte CEM, nous retiendrons uniquement le temps t; qui 
correspond à la plus importante variation du courant lors de l'ouverture. t; est de 
l'ordre de la microseconde. 


L2.4.3 Le MCT 

C'est un thyristor commandé à l'ouverture (figure L5) disposant d'une grille de 
type MOS (MOS Controlled Thyristor). Le contrôle des commutations peut être 
assuré par l'application d'une tension positive ou négative sur la gáchette. Son 
utilisation est ciblée vers les convertisseurs à résonance. 


ir 


K 


4 x A =À 


: Vp 
i 


Figure L5 : MCT 


Le tableau 1.3 présente les caractéristiques de quelques MCT, (tests réalisés avec 
une charge inductive de 200 uH). 

A noter que ` 

- testle temps de monté du courant à la fermeture ` 

- trest le temps de descente du courant à l'ouverture. 
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di/dt à la | di/dt à di/dt 
V (V) Ic (A) t (us) tr (us) |fermeture [l'ouverture] (Max) 
(A/us) c s) (A/us) 


MINEM 


Tableau 1.3 : Caractéristiques sélectionnées de quelques MCT 


Temps de commutation typique que nous retiendrons : une à quelques 
microsecondes. 


L2.4.4 L'IGBT 


C'est un transistor bipolaire (figure I.6) commandé par une entrée MOS (Insulated 
Gate Bipolar Transistor). Il est trés utilisé de nos jours. Ses performances ne cessent 
de progresser et les temps de commutation se rapprochent de ceux du transistor 
MOS de puissance, méme si un courant de queue persiste en fin de commutation. 


"d | 4 
K K 


Le graphique I.1, construit à partir des recueils de caractéristiques donnés par les 
fabriquants, présente, les temps de commutation de différents IGBT. 


[700 ns 
Otr lu 
1600 ns 
[500 ns. 
[400 ns. 
1300 ns. 
[200 ns. 
100 ns 
ns 
ж « « « « « « ` €£ « ж = = = « « LL «4 4 4 4 4 «4 «4 4 « «4 «4 =: 
$ 2 CSS 8 8 8 & A 8 SZ SS 08 3. 89 SS SS 2 BZR кк к 
8 8 8 8 à S E B B Z Z EOS Z Z Z OR Z Z è B Z B Z Z 38 3 8 р 
8 8 8 8 8 8 8 ES 8 8 S & 8 & S & 8 8 8 8 8 & 8 8 8 8 à 


Graphique 1.1: t, et tp de différents IGBT 
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A l'examen du graphique, on voit que les temps t, et t; définis sur la figure 17 
peuvent différer fortement d'un composant à l'autre. 


90 90 


10 10 


t, te 
Figure 1.7 : Ondes types de l'IGBT 


Nous retiendrons : 
- Temps de montée t, : 30 ns à 200 ns. 
- Temps de descente {;: 50 ns à 200 ns. 
(Il s'agit des temps de montée et de descente et non pas des temps de commutation). 


1.2.4.5 Le transistor MOS de puissance 


C'est le composant le plus rapide utilisé en électronique de puissance (figure 1.8), 
cependant, la valeur du courant qu'il peut commuter limite son utilisation aux petites 
et moyennes puissances. Les familles usuelles présentent un compromis classique 
courant /tension, par exemple : (50 A/500 V), (110 A/200 V), (180 A/100 V). 


ip ip i 
J D Cop H D 
Е — Vos G Vos 
= Cps 
3 | 
Cos s 
Figure 1.8 : Transistor MOS de puissance 


Les commutations sont caractérisées uniquement par les temps de montée (t,) et de 
descente (tr) du courant. On devra toutefois prêter attention aux éléments parasites 
des MOS de puissance, à savoir les capacités C; = Сор + Cas et Cps, ainsi qu'à la 
diode de structure en antiparalléle qui peut rendre le MOS bidirectionnel en courant. 

Le graphique 1.2 suivant, construit à partir de recueils de caractéristiques des 
fabricants, indique les temps de commutation de différents types de MOS de 
puissance. 
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350 ns 4 


| 
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Graphique 1.2 : t, et t; de différents MOS de puissance 


On constate à nouveau une disparité importante entre composants. Nous 
retiendrons que les temps de commutation des MOS de puissance rapides se situent 
globalement en dessous de 100 ns, et que les plus rapides d'entre eux peuvent 
commuter en moins de 40 ns. 

L'ordre de grandeur retenue sera donc de 20 à 100 ns. 


1.2.4.6 Le transistor bipolaire de puissance 


Le transistor bipolaire de puissance est ici utilisé en régime bloqué/saturé pour 
assurer les commutations. 

Il concerne principalement les petites et moyennes puissances (commande de 
moteurs, alimentation à découpage faible coût). 

Il convient ici de citer le montage « darlington » qui permet de réduire le courant 
de commande, mais qui n'est guère plus usité de nos jours. 

Des modules intègrent directement ce montage sur le silicium, mais leur diffusion 
commerciale diminue au regard de l'intérêt présenté par les IGBT. 

Le tableau L4 suivant donne les temps de commutation de différents transistors 
bipolaires. 
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Tableau 1.4 : Caractéristiques de différents transistors bipolaires 


Légende pour les tests sur charge résistive : 
ton : temps entre début de la commande et fin de la fermeture du composant ; 
- (В) : temps de descente du courant. 
Nous proposons de retenir un temps de commutation «typique» de 100 à 
quelques centaines de nanosecondes. 


1.2.4.7 Le triac 


Le triac est un composant bidirectionnel en tension et en courant (figure I.9). Il est 
commandé à la fermeture, et son ouverture est naturelle. Le triac est l'équivalent de 
deux thyristors montés en antiparalléle. Sa commande est toutefois dissymétrique, et 
quatre modes de fonctionnement en résultent. 


Аз 


af 7 
Figure 1.9 : Triac 
Les principales caractéristiques (tableau I.5), que nous retiendrons de ces 
composants sont : 
-  ledi/dt à la fermeture ` 
- le dv/dt critique à l'ouverture pour un dig/dt de commande donné ; 
- le dv/dt à l'état bloqué. 
Nous retiendrons que : 
-  lestemps de commutation sont de l'ordre de plusieurs microsecondes ; 
- les dv/dt critiques sont assez faibles (10 V/us) ce qui limite l’utilisation des 
triacs à des puissances moyennes et aux basses fréquences. 
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— 


Ig = 500mA 
Ig = 500mA 
= SUA 
Ig = 500mA 
LAURE 
Ig = 200mA 
i 
á 


5 
2 5 dig/dt = 0,2 A/us 
Ig = 200mA 
5 25 dig/dt = 0,2 A/us 
Ig = 200mA 


Tableau I.5 : Caractéristiques de commutation de triacs de petite puissance 


I.2.5 Comparatif et bilan 


Notre objectif est d'estimer les durées possibles de croissance (tj) et de 
décroissance (ts) du courant principal des composants. Pour cela, nous nous limitons 
au panel retenu précédemment, pour lequel les tensions et courants maximum sont 
comparables. A noter que t, et t; sont estimés sur une plage de niveau du courant 
(10 % - 90 95) ce qui exclut l'intégration des données correspondant aux « courants 
de queue ». Le tableau L6 indique les plages de valeurs plausibles de t, et t; pour les 
échantillons des diverses familles sélectionnées. 


(*) valeur pour des composants de moindre qualité, ou de technologie plus ancienne. 


Tableau 1.6 : Comparatif des temps t, et t; des différents semi-conducteurs 


En observant le tableau récapitulatif, il ressort que deux groupements essentiels 
sont envisageables au regard des temps de commutation définis. 
1* groupement : « thyristor + GTO + triac + transistor bipolaire » ; 
- 2 groupement : « IGBT + transistor MOS de puissance ». 
Par recouvrement approximatif des fourchettes de valeurs, nous pourrions retenir 
pour le premier groupement les chiffres moyens suivants : t, — 0,2 us et te= 0,7 us et 
pour le deuxième groupement : t, # t; # 0,1 us. 
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A ces chiffres, on pourrait faire correspondre des fréquences « équivalentes » 
définies en CEM pour une transition supposée unique par l'expression: 


Fa = ea ce qui conduirait aux indications de valeurs suivantes : 
A D 
(r.f) 


- premier groupement 450 kHz P 1,6 MHz 
- deuxième groupement Fe # 30 MHz 


1.3 Formes d'ondes des courants et modèles simplifiés correspondants 


Suivant la nature des sources et le type de conversion, les formes d'ondes des 
courants seront trés variées. Les éléments passifs, mais également les impédances 
parasites vont influencer fortement les allures. Il est évident que l'on ne peut pas 
traiter le probléme du rayonnement de maniére exhaustive. Sachant que les 
commutations les plus rapides seront les plus génantes au regard des perturbations 
conduites et rayonnées, nous proposons alors de traiter plus spécifiquement le cas 
des cartes de conversion à découpage. 

L'onde de base du courant (sans tenir compte de l'environnement, et en 
privilégiant les aspects haute fréquence) est un trapèze asymétrique ou symétrique, 
ou éventuellement une onde rectangulaire. 

Les ondes particuliéres, isolées, ou superposées à cette onde de base, sont les 
impulsions et oscillations inévitables liées aux inductances et aux capacités 
parasites, des pistes et des SC. 

Avant de passer à la détermination dans le domaine fréquentiel de ces ondes, nous 
jugeons utile de présenter deux montages simples susceptibles de générer ces ondes 
particuliéres. 


1.3.1 Impulsions très brèves de courant d'un pont à sortie 3 niveaux 

Voici le schéma de principe d'un pont onduleur (figure 1.10), fonctionnant à 
20 kHz, équipé de transistors MOS de puissance dont la capacité « drain-source » 
est de l'ordre de 400 pF. La charge prévue est résistive. 

L'onde de sortie est une onde traditionnelle « 120? ». 

La commande retenue par le concepteur se veut la plus économique. Les MOS 
sont contrôlés par paires (une paire = transistor du haut d'un bras et transistor du bas 
de l'autre bras), avec signal de grille appliqué durant 1/3 de période et sans décalage 
du deuxiéme bras sur le premier. On sait que ce choix ne permet pas de forcer la 
tension de sortie à zéro si la charge n'est pas strictement résistive. Ce sera le cas ici, 
puisque nous tenons compte de la présence d'une inductance de connexion en sortie. 

L'apparition de signaux perturbateurs pouvant être à l’origine de l'apparition des 
champs proches est dans ces conditions nettement favorisée. 
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MOS 
observé 
+ 
Source 
d'entrée <+) 


Diode 
73 associée 


Inductance de 
connexion 
Ра р оси 


Charge 
en sortie 


Figure 1.10 : Pont onduleur à transistors MOS de puissance 


Les ondes essentielles relatives aux éléments signalés sur le schéma sont obtenues 
par simulation avec le logiciel EDWIN 32" (figure 1.11). On distingue : 

- des oscillations de tension liées à Ja nature malencontreusement inductive de 
la charge et à l'existence des capacités « drain-source » des transistors MOS. Elles 
apparaissent à l'ouverture, aux bornes des MOS et en sortie du montage. Elles se 
positionnent à mi-amplitude, quand nul semi-conducteur n'est passant. Ces 
oscillations seront génératrices de champ électrique. 

- Des impulsions de courant dans les diodes, à croissance très rapide, de durée 
dépendant de l'importance du stockage inductif. (On notera que durant ces 
impulsions, la tension de sortie change de polarité et que le potentiel du drain du 
transistor MOS rejoint la tension source. Ces niveaux seraient annulés avec une 
commande de bras décalée qui ferait apparaitre, comme il est bien connu des 
électroniciens de puissance, la conduction mixte « transistor - diode », inexistante 
dans le cas présent.) 
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(Reprise par les 
diodes) 
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Figure 1.11 : Ondes essentielles du pont 


Un zoom de l'encadré 1 est présenté sur la figure I.12 dans lequel on notera les 
transitions rapides des courants. 
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Figure 1.12 : Zoom de l'encadré 1 


On remarque ensuite dans l'encadré 2 de la figure 1.12 (détail donné figure 1.13), 
que des impulsions de courant élevées et beaucoup plus rapides encore, apparaissent 
à l'instant de fermeture des transistors d'une paire. Cela est dû à l'application 
brusque qui en résulte de la tension continue d'alimentation sur les jonctions drain- 
Source des autres transistors. En effet, dans l'intervalle précédent de non-conduction 
globale, les capacités parasites drain-source des quatre transistors s'étaient chargées 
à hauteur de 50 96 de cette tension d'alimentation. Il y a donc doublement brutal de 
la tension sur ces éléments capacitifs, sans cheminement pour le courant par 
l'inductance en sortie. 
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Figure 1.13 : Zoom de l'encadré 2 
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Ces dernières impulsions de courant, de niveau double dans la boucle 
d’alimentation, sont de durée très courte, car voisine de 3 à 4 ns. Elles peuvent donc 
générer un champ magnétique d’amplitude significative dont la fréquence 
équivalente peut atteindre 100 MHz. 


1.3.2 Oscillations de courant en amont ou aval d'un hacheur rapide 

La figure I.14 représente le schéma d'un hacheur 100 V/100 kHz. La charge est 
une résistance de 20 О, à laquelle il convient d'adjoindre une inductance de ligne de 
2 uH. L'interrupteur comprend un transistor MOS de puissance de capacité drain- 
source relativement élevée (400 pF) et une diode rapide en anti-paralléle. 


Figure 1.14 : Hacheur rapide 


Des oscillations vont naturellement apparaître, aussi bien en tension, qu’en 
courant. La figure I.15 représente les ondes temporelles qui en témoignent. 

Il faut alors constater que l'oscillation de courant, dans la charge, dont la 
fréquence est voisine de 700 kHz et dont l'amplitude atteint 4 A, sera génératrice 
d'un champ magnétique perturbateur émis par la boucle indésirable que constitue la 
liaison au récepteur. 
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Figure 1.15 : Ondes temporelles du hacheur 
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Un zoom de l'encadré de la figure I.15, est présenté figure 1.16. On peut ainsi 
observer que la diode est sans effet sur l'oscillation de tension dés lors que 
l'amplitude de celle-ci est inférieure à la tension source. 
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Figure 1.16 : Zoom de l'encadré de la figure 1.15 
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1.3.3 Ondes de base dans les domaines temporel et fréquentiel 
Nous allons examiner successivement les différents cas précités, avec pour 
objectif la détermination de leur spectre fréquentiel, lequel sera utilisé 
ultérieurement pour le calcul du champ magnétique émis par des boucles de câblage. 
Les ondes de bases vont être présentées dans le domaine temporel, leur 
décomposition spectrale est obtenue par le calcul des coefficients de Fourier. 


Rappel sur les coefficients de Fourier 
Soit x(t) une fonction périodique de période T. Cette fonction peut s'écrire sous la 
forme d'une somme de fonctions harmoniques : 


e 7 
x(t) = ag DE poi b, db mj (1.1) 


n-i 


avec: 
1 
=— (1.2 
N T [rod ) 
a, [xo * cos it. t (1.3) 
T T 
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b, -2 [xtsin Z, t (14) 
Т; Т 


En réécrivant cette égalité sous forme exponentielle complexe, on obtient la 
décomposition en série de Fourier de la fonction x(t) : 


= ju 22 í 

x)= Yore Т (L5) 
Deen 

avec 
1 - jn, 16 

c, => joe T dt (1.6) 
dëi 

су ear c, _4 - Jb, pour n >0, с, 24. +), pour n < 0. 


Les coefficients c, sont les coefficients de Fourier, mais seuls les c, pour n > 0 ont 
une signification physique. On peut en effet exprimer x(t) comme suit : 


X(t) =c, + y2 


c оз ЕГУ Ze.) (Zc, : argument du complexe c, ). (17) 
n T n 


1.3.3.1 Onde rectangulaire 


io 


tm 
T 
Figure 1.17 : Onde rectangulaire de rapport cyclique œ = ty / T 
EJCD 012 : Onde rectangulaire pour œ = 0,04 (figure 1.17а) 


Pour une onde rectangulaire (figure 1.17), les coefficients de Fourier sont : 


с, = А KA (1.8) 
пл Т 
Les amplitudes des harmoniques sont donc : Н, = 2lc, | (figure 1.18). 
„ = 2A -isin NM ii (1.9) 
nT T 
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Figure 1.18 : Décomposition spectrale de l'onde rectangulaire (rapport cyclique œ = 0,04) 


L'enveloppe des amplitudes des harmoniques décroit suivant une pente de 
moins 20dB/décade à partir du rang T/(ty .n). 

Si le rapport cyclique augmente, l'excursion du piedestal (niveau de départ 
observable sur la figure 1.18) va se réduire et les perturbations liées aux premiers 
rangs des harmoniques seront minimisées. 

Cela est illustré sur le CD d'accompagnement : 


EXCD 014-015 : Représentations temporelle et fréquentielle de l'onde rectangulaire pour œ = 0,18 (figures 
L18aetL18b) 


CD 016-017 : Représentations temporelle et fréquentielle de l'onde rectangulaire pour œ = 0,9 (figures 
L18c et 1.18d) 
1.3.3.2 Onde trapézoïdale symétrique 
ig 
M 


Figure 1.19 : Onde trapézoïdale symétrique 


MCD 018 : Onde trapézoïdale symétrique pour œ = 0,15 et t, = 500 ns (figure I. 19a) 
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Pour un trapèze régulier (figure I.19) les coefficients de Fourier sont : 


D 2. _ Zait 12) 
AT dk T-t, IE л чә), тиын (L10) 


n 


"ray e T 
Les amplitudes des harmoniques sont données par : H, — 2-|C,| (figure 1.20). 


= 2AT = {sin n T-t, MA nT (t, tty) (1.11) 
t, Gn) T P 


n 


20 


Amplitude (dB) ` 


-60 


-70 31 H 
10 10 10 
Fréquence (Hz) C 


Figure 1.20 : Décomposition spectrale de l'onde trapézoïdale symétrique pour œ = 0,15 


L'enveloppe des harmoniques présente trois pentes distinctes : 

- Ја première est nulle pour les harmoniques de rang inférieur à T/(x.(t,+tu)) ; 

= la seconde de - 20 dB/décade pour les harmoniques compris entre 

Ui (tty) et T/(m.t,) ; 

- la troisième de — 40 dB/décade pour les harmoniques supérieurs à T/(m.t,). 

Quand le rapport cyclique augmente, le piedestal disparaît, comme illustré sur le 
CD avec a= 0,45. 


ACD 020-021 : Représentations temporelle et fréquentielle de l'onde trapézoidale pour œ = 0,45 (figures 
1.20a et 1.20b) 
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1.3.3.3 Onde trapézoïdale asymétrique 
iw 


M E 


i t tm tr B 
T 


Figure 1.21 : Onde trapézoidale asymétrique 
ВСР 022 : Onde trapézoidale asymétrique pour œ = 0,2 t, 2 400 ns t; 2 2 ms (figure L21a) 


Les expressions des coefficients de Fourier d'une telle (figure I.21) sont : 


AL А Л EE) 
Cut — wu | 244 SM EES? 
t, t, (2-m- n) T m 
ПЛСР, = BEE f Dm tty) 
oul = d Т T т 


La figure 1.22 donne la décomposition spectrale de ce type d'onde, pour un 


rapport cyclique œ = 0,2 et des temps de transition t, et t& égaux respectivement à 
400 ns et 2 ms. 


20 [34 ï EF TT 


10 


Ф 


ко 
c 


Amplitude (dB) 


ә 
c 


-50 


| 


uU 


Fréquence (Hz) Lio; 


Figure 1.22 : Décomposition spectrale de l'onde trapézoidale asymétrique 


-60 
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L'enveloppe des harmoniques présente trois pentes distinctes dont les points 
d'inflexion sont proches de ceux relatifs au cas de l'onde trapézoidale symétrique. 
Les niveaux sont cependant légérement différents comme indiqué sur la figure. 

En augmentant œ, comme précédemment, le piedestal peut de nouveau disparaître 
et 2 pentes au lieu de 3 subsistent pour l'enveloppe des maximas du spectre. 


АСР 024-025 : Représentations temporelle et fréquentielle de l'onde trapézoidale pour œ = 0,5 (figures 
1.22a et L22b) 


1.3.4 Prise en compte des oscillatio ns de commutation 

Au sein d'un convertisseur à découpage, il va sans dire que des courants de formes 
diverses vont apparaitre, porteurs d'oscillations. Il est bien évidemment difficile de 
vouloir prendre en compte toutes les allures possibles. 

Nous proposons donc de ne considérer que l'allure trapézoidale précédemment 
étudiée, à laquelle nous superposons une oscillation de commutation à l'ouverture 
mais éventuellement également à la fermeture, sachant que ce cas de figure est tout à 
fait envisageable. 

Formulation mathématique 

Les oscillations parasites vont étre modélisées par une sinusoide amortie du type : 


p(t) = Dp-e fi cos(w, «t 9,)- (1.13) 
Représentations temporelles et fréquentielles 


e Cas de l’oscillation isolée 


I (A) 


0 0.5 1 1.5 2 
Temps (s) 0 Eis 


Figure 1.23 : Onde temporelle de l'oscillation isolée, f = 8 MHz 
SCD 027: Oscillation à 8MHz (figure 1.23a), zoom de l'onde de la figure 1.23 
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Amplitude (dB) 


10° 10” 
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Figure 1.24 : Décomposition spectrale de l'oscillation isolée de fréquence 8 MHz 


« Cas de l’oscillation à la fermeture portée par une onde trapézoidale type 


14 
12 
10 | ï 
| | 
ne | | 
2 \ | 
T ob | | 
6 | 1 1 
| | | 
4 | \ \ 
d | \ 
% 0.5 j 


1 
Temps (s) 


Figure 1.25 : Onde temporelle de l'oscillation à la fermeture de fréquence f = 8 MHz 
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Figure : 1.26 Décomposition spectrale de l'oscillation à la fermeture de fréquence f = 8 MHz 


e Cas de l’oscillation à l'ouverture portée par une onde trapézoidale type 


1 
Temps ($) 
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Figure 1.27 : Onde temporelle de l'oscillation à l'ouverture de fréquence f = 8 MHz 
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Figure 1.28 : Décomposition spectrale de l'oscillation à l'ouverture de fréquence f = 8 MHz 


Commentaire 


П apparait clairement que la présence d'une oscillation parasite modifie 
l'enveloppe du spectre du courant de maniére significative aux fréquences les plus 
élevées. Le phénoméne sera amplifié par dérivation de la fonction comme nous le 
prendrons en compte ultérieurement. 


Conclusion du chapitre 


Nous venons de voir briévement ce qu'est la fonction de commutation en 
électronique de puissance, et quelles sont les structures de conversion classiques qui 
l'utilisent. Nous avons également présenté une étude statistique des temps de 
transition d'un échantillonnage de semi-conducteurs associés à ces structures. 

Deux groupements de semi-conducteurs peuvent étre retenus. Les plus rapides 
présentent des temps de commutation conduisant à des fréquences équivalentes 
pouvant atteindre 3 MHz. On peut donc s'attendre à un rayonnement 
électromagnétique important au voisinage des circuits associés. 

De maniére à aborder ce probléme, nous avons sélectionné des ondes typiques de 
courant, construites à partir des considérations précédentes sur les temps de 
transition des semi-conducteurs, mais également sur la présence possible 
d'oscillations parasites liées aux imperfections des cáblages. 
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CHAPITRE II 


PRÉSENTATION SIMPLIFIÉE 
DU 
RAYONNEMENT 
DES 
CARTES DE PUISSANCE 


Les convertisseurs réalisés sur cartes sont de petites et moyennes puissances, mais 
ils constituent le plus grand nombre. Nous allons donc cibler essentiellement notre 
présentation sur ces structures coplanaires 2D, sachant que les systémes de grande 
puissance nécessairement tri-dimensionnels doivent faire l’objet d’études 
spécifiques individualisées. 


11.1 flot rayonnant, boucle « topologique » et boucle « équivalente » 


11.1.1 « flot rayonnant » 

Les mesures du champ magnétique au-dessus de cartes électroniques, que l'on 
peut effectuer à l'aide de sondes de champ proche, permettent en général de repérer 
des « flots rayonnants » localisés. 

Pour une carte simple d'alimentation à découpage par exemple, on pourrait repérer 
deux ou trois « 1015 » localisés prépondérants dont on peut considérer que l'action 
conjuguée détermine l'essentiel du champ magnétique résultant de la carte. Dans ces 
conditions, nous allons dans un premier temps étudier tous les aspects du 
rayonnement de boucles particuliéres à géométrie simplifiée et conjuguer leurs effets 
dans un deuxième temps. 


11.1.2 « Boucle topologique » 

La « boucle topologique » voit son nom inspiré par la notion de « topologie » qui 
est une forme de représentation du schéma structurel limité aux semi-conducteurs et 
pistes parcourus par des courants. La « boucle topologique » est donc une partie de 
la structure à un instant donné. Nous allons expliciter cette notion sur quelques 
exemples. 

La présentation qui va étre faite se veut didactique et se base sur des schémas 
simplifiés, auxquels peuvent correspondre des implantations géométriques voisines. 

Notre objectif à terme est d'obtenir le rayonnement de « boucles équivalentes » 
aux flots précités, dans l'espace environnant (3D). 

Hacheur abaisseur 

Le schéma retenu est donné par la figure II.1. 

Dans cet exemple, deux topologies essentielles peuvent étre sélectionnées, 
auxquelles correspondent les deux « boucles topologiques » suivantes : 

- Bu [С, T,Z]; 

- Bu:[D,2] 


En implantation réelle, la boucle By peut être constituée simplement par les pistes 


de sortie vers les bornes de connexion de la charge, souvent reliée par un câble 
bifilaire ne présentant pas de surface de boucle. 
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Figure II.1 : « Boucle topologique » d'un hacheur abaisseur 


La présence d'un CALC associé au transistor T peut évidemment ajouter d'autres 
« boucles topologiques ». 


Pont monophasé 
Le schéma de principe d'un onduleur de tension est donné figure II.2. 


Figure II.2 : « Boucle topologique » d'un pont monophasé 


L'implantation réelle pourrait étre de géométrie voisine de celle du schéma, ou 
trés différente suivant, par exemple la disposition du (des) refroidisseur(s), la 
présence de potentiels communs (exemple celui des émetteurs des semi-conducteurs 
du bas) etc. 

Notre présentation ne pouvant étre exhaustive, nous proposons d'observer les 
topologies sur le schéma de principe. 

Les « boucles topologiques » Bu et B, ont pour particularité d’être : 

- de méme surface ; 

-  symétriques ; 

- parcourues par des courants de mêmes sens. 

Elles présentent en outre une partie commune incluant le condensateur haute 
fréquence C. On notera que la charge supposée ramenée comme précédemment par 
un cáble bifilaire, n'introduit pas de boucle extérieure. 
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IL.1.3 e Boucle équivalente » 


Les logiciels de routage proposent en général des parcours de pistes tenant compte 
des encombrements des composants, des densités de courant et de règles essentielles 
pour limiter les effets de diaphonie et de couplage par impédance commune. Les 
parcours des pistes sont parallèles ou perpendiculaires. Les « boucles topologiques » 
qui peuvent se dégager sont donc essentiellement des formes géométriques 
construites à partir de segments perpendiculaires. 

Sur une carte réelle, les «boucles topologiques » correspondent à des semi- 
conducteurs, des composants passifs, et des éléments de piste déterminés. Dans leur 
environnement, il existe également des composants à faible courant destinés par 
exemple à la commande des interrupteurs, à la régulation, au découplage etc. 
Sachant que les «boucles topologiques» constituent l'essentiel du champ 
magnétique rayonné à proximité, on peut alors cerner « géographiquement » un îlot 
correspondant. 

Une carte standard d’alimentation à découpage présentera plusieurs « îlots », donc 
sera schématisée par autant de « boucles équivalentes » (B.) de forme simples. 

Le rayonnement global de la carte observée au travers de celui des « boucles 
équivalentes » sera estimé au sens de la CEM, c'est-à-dire dans les conditions les 
plus défavorables. 


II.2 Formulation mathématique et représentation spatiale du champ 
magnétique 
La plupart des ouvrages présentent essentiellement la boucle élémentaire 
(circulaire infinitésimale). Il convient de s'interroger sur plusieurs points, 
concernant la comparaison possible de cette boucle élémentaire avec la « boucle 
équivalente » réelle type de surface finie, de forme rectangulaire ou carrée. 
Voici les points d'interrogation qu'il faut retenir : 
- quel est le domaine de validité en fonction de l'éloignement, de l équivalence 
possible entre la boucle élémentaire et la boucle rectangulaire ? 
- quelle est la validité d'une démarche de simplification d'écriture 
mathématique basée sur le paramètre Br (produit : nombre d'onde x distance) ? 
- quelle est l'influence de la dimension des conducteurs constituant la boucle ? 
- comment prendre en compte l'effet d'un plan de masse éventuel situé 
immédiatement sous le circuit ? 


IL2.1 Equations électromagnétiques 
Expressions locales des équations de Maxwell : 


УХЕ РЕ 
ot Wf 


VxH= —— (1.2) 
x J+ E 
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V.D-p (11.3) 
V-B=0 (11.4) 
Précisions sur les opérateurs et grandeurs considérées : 

-  V-: calcul de divergence ; 

- Ух: calcul de rotationnel ; 


: vecteur champ électrique (V/m) ; 

: vecteur champ magnétique (A/m) ; 

: vecteur d'induction magnétique (Vs/m2 ou T) ; 

: vecteur champ de déplacement (As/m? ou C/m?) ; 
: vecteur densité de courant (A/m?) ; 

-  p:densité volumique de charge (C/m’). 


SS Do Do 8 


Dans les milieux isotropes homogènes nous retiendrons les relations suivantes 
entres les grandeurs : 


B= uH - uu, H (11.5) 
D=e€É=e,€,Ë (IL6) 
J=oE (11.7) 
avec : 


- £:permittivité absolue, £r : permittivité relative, £0 = 1/(36.7.109)As/Vm ; 

- p: perméabilité absolue, u, : perméabilité relative, ho = 4.x.10" Vs/Am ; 

- ©: conductivité S/m. 

Du fait de V. B = 0, il existe un vecteur potentiel magnétique À tel que : 
B-VxA (11.8) 
Ce vecteur potentiel n'est pas complètement défini par cette relation, il le sera 

ultérieurement par l'expression de sa divergence. 

En introduisant (П.8) dans (II.1) on obtient : 


- дА 
Vx| E+ |=0 ш) 


Comme le rotationnel d’un gradient est nul, cela permet de définir le potentiel 
scalaire électrique V suivant : 


E+-—=-V(V) (11.10) 


V(V) étant le gradient de V. 


Ce potentiel scalaire électrique est alors défini à une constante près. 

A partir de l'équation (2) et en utilisant (8) et (10) on obtient : 

VxVxA=uJ йе A e (1.11) 
d'r ot 
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En utilisant la relation Vx Vx À = V(V.A) - V/A, (11) devient : 


= дА > 7. 

V?A-pe—— = -uJ + V| V.A + u£ — (1.12) 
H Бүл H | H ot ) 
En utilisant (10) dans l'équation (3) il vient : 
dE (11.13) 
дї E 
Sachant que V-V(V) = V?V , on peut écrire (13) sous la forme : 
dus 
Vj ef. pod Т a (11.14) 
971 € at ot 


Pour définir complètement le vecteur potentiel magnétique et le potentiel scalaire 
électrique on impose la condition de Lorentz : 


= ov 
V.A+ u£— = 0 (1.15) 
а 
Il s'en suit que (12) et (14) deviennent : 
Е д?А 2 
V'À-uc2 A Al (116) 
H 35. H 
= av p 
Vy — ue = = Er (1.17) 
py E 


Ce sont les deux équations d'onde. 
Les solutions particulières de ces équations sont : 


AM t) = H Í es (11.18) 
"Te F 
aah 2525 _ * (11.19) 
"cr 1] = td | 


M : point d’observation 
R : distance entre le point d'observation M et l'élément dv 
Ces solutions vérifient les conditions aux limites suivantes : 


lim À = 0 €t іту =0 


re ғ 

А(М ,t) et V(M,t) sont des potentiels retardés. C'est-à-dire que leurs valeurs à 
l'instant (t) dépendent du courant ou de la charge à l'instant (t-r/u) correspondant au 
retard dû à la propagation à la vitesse u des effets électromagnétiques. 

Nous allons maintenant uniquement nous intéresser au vecteur potentiel 
magnétique. 

Tout signal périodique pouvant étre décomposé sous la forme d'une somme de 
fonctions sinusoidales, nous pouvons introduire la notation complexe dans l'écriture 
des formules. 
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Soit G(r) un vecteur. En régime sinusoidal permanent et en coordonnées 


cartésiennes il peut s'exprimer sous la forme : 
G-2G,X4G, 34 G.Z 
avec: 
G. = Ср, cos(wt + Фф, ) 
G, = Guy cos(wt +, ) 
G, =Gy, cos(wt+ @, ) 
En utilisant l'écriture complexe on obtient : 
G, = Se(G,, ее!) есе") avec: G, = Gye!” 


G, = Re(G,,e"e” )=RelG,e™ ) avèc : G, = Gy e 
G, = e(G,, e^" e^" )- Zelt e?" ) ауес: G, = Gye!” 


x 


(11.20) 


(11.21) 
(11.22) 
(1.23) 


(11.24) 
(1.25) 


(11.26) 


G,; G,» О, sont des grandeurs complexes indépendantes du temps (parfois 


appelées phaseurs) dont le module représente l'amplitude maximale de G(r) suivant 


chaque direction. 


A partir de ces grandeurs, on peut définir le vecteurG (parfois appelé vecteur- 


phaseur) : 
G-G,X4G, 3 *G.Z 


La relation qui relie G et G(r) est donc : Gu) = Re(Ge™ ) 
Remarque : 


200 «ence 


(1.27) 


(11.28) 


(11.29) 


Le vecteur-phaseur sera souvent utilisé dans la suite de cette présentation car il 
permet de déterminer l'amplitude maximale du champ tout en simplifiant les 


écritures. 


A partir de cette écriture on peut exprimer en régime sinusoïdal les expressions 


complexes du courant excitateur et du potentiel vecteur : 


j= Relje оо) 


- J: vecteur densité de courant dans le conducteur ; 
- J: T densité de courant. 


Ja ——) 
АМ) = 7 re и reg 
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11.2.2 Application à la boucle élémentaire 


II.2.2.1 Formulation mathématique 
Le milieu de propagation est l'air, donc р = poet € = £o. 
Hypothéses de calcul : 
- les conducteurs sont supposés filiformes donc : 


ff Jdv = fiai (11.32) 
i 1 


- les dimensions des boucles sont petites devant la longueur d'onde des 

fréquences considérées. Le courant est donc constant tout le long de la boucle. 

Nous considérons une excitation sinusoidale du courant dont nous prenons 
l'expression complexe : 

I I0" (11.33) 

L'expression complexe du potentiel vecteur magnétique est donc : 


i(wt—Br) 
EE L (11.34) 
4л : r 
avec 


B- 27 nombre d'onde 


À : la longueur d'onde associé à la fréquence f 
L'écriture simplifiée (sans faire apparaître la dépendance temporelle) est donc ` 


А= Ё. ku (1135) 

"Arr 

La détermination du champ magnétique se fait à partir de (IL5) et (IL8) : 

e. VXA (11.36) 
Ho 


Pour une boucle élémentaire (figure II.3), les composantes des champs sont en 
coordonnées sphériques. 


nus z+ ES (cos0)e ^? (11.37) 
E 2л f (Br) (Br) 


= p _ ! j sin @)e Jr (11.38) 
Ho SE (Br) ( (Br? H пӨ)е 
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Fig.Il.3 : Système de coordonnées pour la boucle élémentaire 

avec: 

- I: amplitude du courant excitateur ; 

- ` Se: surface de la boucle élémentaire. 

Suivant la valeur du produit fjr on distingue deux régions de l'espace ou l'on peut 
simplifier les écritures du champ magnétique. 

En champ proche : Br << 1, les termes dominants sont ceux de plus haut degré ; de 
plus on peut considérer que e GP? #1. 

Dans ces conditions, les composantes du champ magnétique deviennent : 


H, Be sin@ (IL.39) 
* Amr 

Н, = e cos (11.40) 
2л 


L'amplitude du champ magnétique est donc ` 
_ IS, [3cos(20)+5 

2nr? 8 
Donc, en champ proche, l'amplitude du champ magnétique est indépendante de la 


fréquence. Elle ne dépend que du courant dans la boucle, de la surface de celle-ci et 
du point d'observation. 


La frontière de cette zone de champ proche sera précisée plus loin. 
En champ lointain : Br >> 1, ce sont les termes de plus bas degré qui dominent. 


(11.41) 


3 
H, BP [s fino K (142) 
— 4m | Br 
3 E 
H, = EX 5 j= KE (11.43) 
Л r 


Hoest en (1/Br) et H, est en (1/Br)°. Dans ces conditions on ne considère que Hg 
Dans le cas général : 


3 . 
Bi Rele" )= IS,p [se - Br) cos(wr- Br) sin(wt- Br) EN (11.44) 


л | (ве) (Br) (Br) 


40 


Chapitre II Présentation simplifiée du rayonnement des cartes de puissance 


= B jwt _ 158? cos(wt— Br) _ sin(wt - Br) (11.45) 
EE | (Bf. @ Fr" 

8 DCN | ga SinQwr- Br) [sine | 
H = SE (By (Br) cos(wt — Br) (ey Э 


(Br) (Br) 


П.2.2.2 Représentations spatiales 


L'objectif est de mettre en évidence la distribution du champ en tout point de 
l'espace sur la base de trois descriptions particuliéres : 
(a) Le lieu des points est un hémisphére recouvrant la boucle élémentaire. 
(b) Le lieu est constitué par des plans paralléles, tous perpendiculaires au plan 
de la boucle et réguliérement éloignés de celle-ci. 
(c) Le lieu des points est un quart de plan perpendiculaire à la boucle 
élémentaire et dont l'intersection avec celle-ci est un rayon (localisation en 
2D par r et 0). 
Représentation suivant la description (a) 

La figure П.4 donne l'amplitude maximale et l'orientation en tout point de 
l'hémisphère, à fréquence donnée, du champ magnétique proche. Cette 
représentation permet de souligner que les directions du champ ne sont évidemment 
pas radiales (sauf à 6 = 0) mais résultent d'une composition de vecteurs ( H 5 Н 2d. 


1 élémentaire |. _ _ - N É dn ------- Paramètres: 
| centré : F = 100 kHz 
centrée en MM Y > > < 
99У 4000) XX Witty, 000 S=4cm 
CRD Ce Ei DUE E I 210A 
ч кык зш ee e Tsphére = 1 M 


X (m) 
Elo 
Figure П.4 : Champ magnétique sur un hémisphère centré sur la boucle élémentaire 


SCD 034 : Représentation pour une sphère complète (figure IL4a) 
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Nous allons examiner raaintenant l'influence de l'éloignement (r) sur la variation 
d'amplitude du champ magnétique en considérant celui-ci dans le plan Y-Z. 

Remarque : nous sommes en « champ proche », on pourra vérifier la concordance 
avec les expressions mathématiques sur le rapport 2 entre les amplitudes observées à 
la verticale de la boucle et dans son plan. 


0.12 F 


` D» 


0.08 | Ae 


` 
e 0.06 + ERR 
< >> 
Paramètres : à. 
0.04 F F =1 MHz 
Se = 4сп? = 
0.02 - 1=10А NX 
oF \ 
Е d. 
-0.02 0 002 0.04 0.06 0.08 0.1 0.12 
as Eos 


Figure ILS : Champ magnétique proche de la boucle ( r = 10 cm) 


120r 17 
100! 1772 
TT 
80r — > 
ES 
_ 60! E 
& End 
N 40r Paramètres : `` 
F = 1 MHz 
20+ Se = 4 cm? 
Iz10A 
OF 
-20r 
0 50 100 
T.) Eos 


Figure 11.6 : Champ magnétique éloigné de la boucle (r = 100m) 
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Commentaire : nous sommes cette fois en zone de transition et le rapport précédent 
diminue, devenant inférieur à l'unité. 

NB : les conditions d'affichage sont différentes dans un souci de lisibilité. Il 
convient évidemment de ne pas effectuer de comparaison directe d'amplitude. 
Représentation suivant la description (b) 

La figure IL7 est une illustration de la décroissance rapide de l'amplitude 
maximale du champ magnétique avec l'éloignement dans deux directions 
orthogonales (X et Y) pour différentes hauteurs par rapport au plan de la boucle. On 
observe également les changements d'orientation correspondants. 


BP = se sers pee PT EE 
Dares nel ME By ND o ee АШ 
Paramètres : bee oe sy, ' i : 
0.3 .1 2. ne `> Rer LP er É ' 
F =100kHz ` | 799.7 М СУЗ. 
024 S.-4cm .,.... ES . , 
] 210A it DA US ] ч 
pales sedi . 2% fe ES GE 
E УЫ £043 Rit. n LAN e Der E dert | 
— : ч | Gus ' 
Б OUT = Lean Re ^ n AUS Green : 
ЕЧ oc, Se Se BA Uu а Lh 
> уу. О Néi ; 
Er Жн PME : à , 
' WW KN PE e 
VARS 3 ieu ыш d Dess se niodo et A o. `u ' 
MEE SR euren Adi e 
- +» Ï e see Dee Ee, =| 
0 | T = = нара 
-0 0.2 0.4 0.6 0.8 1 
Ma ш) Elo; 


Figure II.7 : Champ magnétique suivant différents plans parallèles à la boucle 


ECD 038 : Vision périphérique du champ magnétique suivant différents plans parallèles à la boucle 
(figure IL 7a) 


Représentation suivant la description (c) de l'espace 

La boucle étant circulaire, on peut ne retenir que les variations en r et 6, et se 
situer ainsi dans un plan perpendiculaire à 1а boucle passant par son centre. 
MCD 039 : Amplitude du champ Н en fonction de (r et 0) à ф = 0 °, avec pour paramètres : f = 100 kHz, 
Se= 4 cm, I=10 A (figure IL7b) 

Bien évidemment, on vérifie que l’amplitude maximale du champ magnétique 
décroît rapidement avec l’éloignement (la fréquence n’étant pas très élevée). 

D'autre part, l'influence de 9 peut être observée aux limites (très près ou très loin 
de la boucle) faisant ainsi clairement apparaitre des sens de variation opposés. 
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11.2.3 Application à la boucle rectangulaire 


Nous jugeons utile d'intégrer dans ce chapitre le calcul complet des expressions 
du champ dans le cas de la boucle rectangulaire qui donnera lieu ensuite à une 
implantation sur MATLAB™ 

Une telle boucle est représentée figure II.8, avec le système d'axes associé. 


Figure IL8 : Boucle rectangulaire 


П.2.3.1 Formulation mathématique 


Pour effectuer le calcul du champ émis par une boucle rectangulaire nous allons 
décomposer l'équation (11.36) en une somme d'intégrales. Chaque intégrale 
correspond à un coté du rectangle dans ce cas. (Si la forme géométrique 
correspondant à la boucle émettrice est plus complexe mais toujours formée de 
SEMIS il suffira еа GE és par pip 

Wi -iBr _ 
4 а „еме Sek See = di, 


(IL47) 
Le vecteur potentiel s'exprime suivant les coordonnées cartésiennes de la facon 
suivante : 


ARANAN A (11.48) 


donc : dA 2 dx +3 zs SÉ + s dz (11.49) 


On notera que E = 0 car la boucle émettrice est dans le plan z = 0. 


Nous obtenons donc : 
; D = jBr = ) -Ц и ) 
A =#0 f e 


x 


— 4л Va 
-4 (2) 
2 


(11.50) 
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& pu) e PL 
ші} te fe (IL51) 
de Í Von dy+ | T -u dy 
= «|t g La 
avec: к= dan = х) + (ум - yV + (zy - 2? (11.52) 
1 
4 = jB[ ом o o e! 72 21 -jB|Gu Hoy Dea! 2 
ш | pe e (11.53) 
Ages Í r dst J > dx 
л zh [ом o - 9+] 4 [Gu tow e! 2 
d d 
сше" [оома ` чё [мом tz E (11.54) 
À, E 4 | p-dy* J 1 dy 
— л = СИИР Ka Ian os ]2 
Calculons le champ magnétique H : 
x y H 
- 1 - 1 d d d (11.55) 
H=—VxA=—det 
Ho Ho Oxy дум 9м 
A А, À 
dA 
Sad (11.56) 
Hx 9 
HO 0м 
PEE: (1.57) 
— ро у 
d 
г [94у däs (11.58) 
Z 


Dans le calcul des expressions des composantes du champ magnétique, on peut 
remarquer que les opérations de dérivation et d’intégration portent sur des variables 
différentes. Il est donc possible de permuter ces opérations. 


9A 
Exemple sur le calcul de _— : 


ER? 
s Lia. du ten 
^y Ə Hof e (11.59) 
dy dm de (ET at LA dr 
2 
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En permutant les opérations on obtient : 


Ural F La 


dA 2 2 2 
59 Hol 9 |е „ [ 9 |e `” (11.60) 
dy 4 à ; "JS r À 
zm 4т |4 OO You wi Tun 
2 2 L 2 
Es ES 
-jfr -L 
3 [=] 
pee = (11.61) 
"Sai T= 
1 
А L2 2 He 
-jB| (xy +? «Gy =>) 2 
с д р D Хм +97) TUM Y) tzM | (11.62) 


i = 
dzy 1, 
D +)? ум tp 


Décomposons cette expression : 


1 
9 | [оом ум | 
=| e 


1 
=—jBzy [Cu +) + (Ju = у) + zu” | "x 


dy 
1 
= BL Gu eise ум P (1.63) 
€ 
et : 
l 3 
12 2 21 5 == 
d EH +(ум У) +2 | 2 XN. EHNEN 2 (11.64) 
Sr = ZM 
dzy 
donc 
RS — JBL Gu eon e! |2 


-jB 1 (11.65) 
~ 3 
[ом +4) Ум wp 


[ом +44 Ум STE 
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Cette expression peut être mise sous la forme simplifiée suivante : 


Ta = tm ; 63] - jp | (11.66) 
z 2 3 
r r 
al ten 
de même : 
— jB | 
д le | (+) WË - jB 1 (11.67) 
Та i d Vu M 2 3 
ч (=t) Lo 


L'expression de Н, est donc la suivante : 


2 І 
B tds LCR P Gy nre D JB T 1 
Hx ie м е ПОГАШ E 3 H 
2. [ob fO XY in ] [ом 0м nte! P 
E ] 
-jB| Gu + Cu —y ten? |2 j 1 
à Je [0 o nt] jp ч - lay} (бв) 
_ 42 TERMS ` Sec 
EEN [о Pow e! E 


Les calculs pour déterminer H, sont similaires, on obtient donc : 


D 1 
2 H 4342 SE: ` 
-iP ом +04 Y +cu P ip 1 
H =-#zy X Í e” : + y 
; = ал м Я $75 3 
= @ zu [ew х7 +0 Bra [ow хум р 
zu 1 ; 
* Í p Momo jp * = Шур (169) 
2 12,2 
л [ow =x) +( yy +P) +2] [ow Cu TEE 
Pour le calcul de H, nous déterminons T,, : 
D 
Е Li e L p 
x=- —jB| (хм +124 (y - ш 2 
2 2 
T= 9 |e _ 9 le (11.70) 
xl — Se 
Oxy (8) dru 4,2 2 2 ! 
2 [ом Жз) +(ум Y) +2м | 
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1 
; L2 2 2h 75 
д | JB Gat) Og У) + 2 f L L 2| 2 
| ? | =-/8\хм t| (ки 2 Hl +(ум у) + zy x 


Oxy 
(11.71) 


-8| ом D y -» tay dd 
e 


et: 


3 
= D A, SEN з = -hru «tes а) ey D ea! * qun 
M 


donc : 
Ta = (xm SC De pow et |2 x 
(11.73) 
- iB z 1 
pou SCH ten? | E +52 HEITEN +ы?|? 
que l’on peut exprimer sous la forme simplifiée : 
"Jr on E 
isis gas [= à) E 5 T и f (11.74) 
(=+) (=) 
de même : 
- Br, . 
EAR (Fal P n) Ta - jB 1 (175) 
T dr "on aide? € а 27 3 
кы (5) te) 
-jbr o. 
E 2) JBr, E e 
Se kal ka 
- j Br E Ls 
r. = ò |e (2) «bu -®), 7-2) -B _ 1 | (1.77) 


y2 dy E 
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On obtient alors l'expression de H;: 


EI 
Po -jpu HE Hy asw 2 
DEE HAN 
m E 
2 


-jp 1 P 
[ow +2) Hyu -rPe | [a e «Ou „йыл 


E 1 


-jB| Gu BY + Cu У) rz |2 
+ (b= gf enter 


-iB dy 
С? -1H yy =y) E `] [E Жум У) +z wE 
ц d 
ES jj e (ay + Ou FY + wp x 
H 
2 
- jB ые M dx 


3 
-х)? 12) 42, ud 
[eu 0 Ou У "e ] [ом +m -2P + | 


Au 1 
— j B| (x, -x)?+( vy +2)? +2 2/2 
(vu +42) fe D | "Py 
4 
2 


=} 1 
iB E - |dx (11.77) 
[Geo ea] 


[ow х) (ум EE 


w 


Les différentes expressions obtenues peuvent être implantées sous un logiciel 
comme MATLAB™ ce qui permet de connaitre en tout point de l’espace le champ 
émis par la boucle rayonnante rectangulaire. 
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П.2.3.2 Représentations spatiales 

La boucle précédente dite élémentaire donne lieu a une formulation mathématique 
simplifiée n'autorisant pas l'examen de l'influence de la proximité du conducteur la 
constituant, puisque seul le produit S,xI (surfacexcourant maximal) est considéré. 
En revanche, dans le cas maintenant de la boucle rectangulaire, le courant Ï est 
parfaitement localisé dans le conducteur (supposé filiforme) décrivant un contour 
déterminé. La figure II.9 permet d'illustrer ce préambule. En effet, si l'on recherche 
le champ H (amplitude maximale et orientation) dans un plan paralléle à la boucle 
au-dessus de celle-ci, on voit se « dessiner » le profil du contour. 


0.035 + 
е Aide / me | 
"E NS Al PALLA i Paramètres: + 
0.025 + AN | | Fra “o į F=100kHz | 
"a ee E Ee uen | = | 


Í x 
fei / j А S 
E 0.02-] 
N 
0.015 - 
0.01 + Wii 
í 0 
KEN Y (m) 
-0.05 | 
-0.1 
0.05 
X (m) Elo 


Figure 11.9 : Champ H dans un plan au-dessus et parallèle à la boucle (Z = 2,5 cm) 


Remarques : 


On peut effectuer des rotations en 3D comme représenté sur le CD : 
FCD 041-042 : Variantes 3D de représentation du champ H dans un plan au-dessus et parallèle à la 
boucle pour z = 2,5 cm (figures П.9а et IL9b) 


L'incidence du rapprochement est évidemment forte : 
WCD 043 : Champ Н dans un plan au-dessus et parallèle à la boucle pour z = 5 mm (figure II.9c) 


Il est intéressant de représenter les courbes de méme niveau d'amplitude du champ 
magnétique, dans le plan de la boucle (figure II.10). 

L'influence de la proximité des conducteurs est d'autant plus faible, que l'on se 
place dans un plan de plus en plus éloigné de la boucle. On pourra observer cela sur 
le CD. 

П est également intéressant de représenter les courbes de méme niveau 
d'amplitude du champ magnétique dans un plan perpendiculaire c'est-à-dire dans le 
plan YZ, en coupe suivant OX (figure II.11). 

L'influence des conducteurs s'atténue cette fois avec l'éloignement suivant OX : 
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0.05 


Y (m) 


0.005 0 0.05 
ш) Es 


Figure II.10 : Courbes d'égale amplitude du champ Н dans le plan (XY) 


ECD 045-046 : Courbes d'égale amplitude de H dans un plan situé à z = 2 et 6 cm (figures H. 10a et b) 
Di 
0.05 + 


Z (m) 


0.05 
Y (m 
» Bos 


Figure П.11 : Courbes d'égale amplitude du champ H dans le plan (YZ) à X = 0 
ECD 048-050 : Courbes d'égale amplitude de H dans un plan YZ situé à x = 2,5 cm; 3,75 cm et 5 ст 


(figures П.11а, IL 11b et IL11c) 
On notera que le champ est bien plus important à proximité des conducteurs, en 


particulier à l'intérieur de la boucle. 
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Si l'on observe la coupe (figure II.11), sachant que le courant est entrant à gauche 
et sortant à droite (ou l'inverse) les champs autour des fils s'ajoutent dans la boucle 
et se retranchent pour un point choisi à l'extérieur. 

Les figures П.12 et П.12а qui montrent le champ trés près de la spire ou du 
conducteur permettent de confirmer directement cette remarque. 


0.05r 
Paramétres : 
i F=100kHz | 
SE ` ï WM =12=5cm | 
dur En epit is i I=10A 
CR C RR RO р с 
m 7 "NOME 03 bt 04 4 4 € < 
E OF VE SE DENT Pr Ў 
м РТС aa 
PRET E E AUX St 
tech v D D D П 1 1 9 ` ` Ki 
-0.05 — 1 1 1 1 1 1 
-0.06 -0.04 -0.02 0 0.02 0.04 0.06 
Y (m) 
Eos; 


Figure П.12 : Champ H prés de la spire 
ECD 052 : Champ H très près du conducteur (figure П.12а) 


` 


Il est possible d'établir une représentation en trois dimensions à @ donné de 
l'amplitude maximale du champ magnétique en fonction de r et Ө (figure 11.13). 

On voit ainsi pertinemment que l'accroissement marqué d'amplitude est le témoin 
de l'influence du conducteur proche. 

La figure IL14 donne l'amplitude du potentiel vecteur A, dans les mémes 
conditions d'observation que pour la figure IL11. Cette nouvelle représentation 
souligne à nouveau de manière flagrante l'incidence de la proximité du conducteur. 

Il convient de rappeler que A est lié au champ par la relation : 


УХА 


Ho 
Remarque sur l’influence éventuelle de la fréquence 
Avec les équations simplifiées de la boucle élémentaire sous l'hypothèse du 
champ proche, nous avons montré aisément l'indépendance de l'amplitude du 
champ de la fréquence. En ce qui concerne la boucle rectangulaire, la mise en 
évidence du phénomène est moins directe, c’est pourquoi nous présentons des tracés 
comparés (figure П.15) pour la boucle élémentaire et la boucle rectangulaire. 


H = 
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H-f(R,0) 


Paramétres : 

F=100kHz : 

5cmx5cm i 
I=I0A | 


Eos 
Figure 11.13 : Amplitude maximale du champ magnétique en fonction de r et 0 (ÈP = 0°) 
SCD 054-055 : Variantes de la figure II.13 (figures II.13a et II. 13b) 
x 10° 
1 
0.8 Paramètres : 
F = 100 kHz 
A 0.6 5 cm x 5 cm 
5 1= 10А 
БУ 
< 0.4 ——] 
0.2 
04 
0.04 2 
0.02 
-0.02 » Я à 0 0.04 
Zu. Oe wp тег T 
ose 


Figure II. 14 : Amplitude du potentiel vecteur À 
ECD 057-058 : Variantes en 3D et 2D de la figure II.14 (figures II. 14a et П.14Ь) 
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Conditions retenues pour effecter les tracés comparés de la figure II.15 : 

- choix du méme point d'observation (X = 0, Y 24cm, Z = 2,5 ст); 

- on évite de situer le point d'observation à proximité immédiate du 
conducteur (cas de la boucle rectangulaire) ; 


- оп fixe le produit S,xI à la méme valeur pour les deux boucles. 


1.05 T - — — = F tt 
` | "Ï — Boucle carrée 
+ Boucle élémentaire |, 


ke 

2 Paramétres : i 

30.95 1стх1ст ecol д 
& 2 1=10А | 

3 + x=0 y=4cm i 

2 i z-25cm | | 
E 0.9 T 
a ] 
= 

Ù 


0.85 


0. 8 A ere Echt d bo EURE ЭРШИ WA RAUS ЖНС titi) 
10° 10° 107 10° 10° 
Fréquence (Hz) Loo 


Figure П.15 : Evolution de l'amplitude du champ Н en fonction de la fréquence 


On remarquera que le changement du champ magnétique ne s'effectue qu'à partir 
d'une valeur parfaitement identifiable pour laquelle on quitterait les hypothèses du 
champ proche. Cela souligne à nouveau l'indépendance de l'amplitude du champ 
d'avec la fréquence (dans la limite d'un éloignement modéré bien sûr). 


IL.2.4 Domaine de validité de l'équi valence « boucle réelle-boucle élémentaire » 

Nous allons examiner la possibilité d'assimiler, au regard du champ magnétique 
émis, une boucle de petites dimensions à une boucle élémentaire. 

Cette équivalence se fait sur l’hypothèse de la conservation de la valeur du produit 
IxS (courant x surface de la boucle). Il n'est bien sür pas concevable d'effectuer cet 
examen pour une boucle réelle de forme trop allongée (pistes paralléles par 
exemple). De plus, les cas concrets d'implantation conduisent statistiquement dans 
le cas du rectangle à un rapport longueur/largeur inférieur à 4. 

Pour des formes plus complexes, le rectangle à coin tronqué (présentant une 
enclave) de la figure II.16 est un cas représentatif que nous pouvons retenir. 

Pour dégager les conditions d'équivalence, il convient de procéder à plusieurs 
études comparatives. 
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Figure П.16 : Piste rectangulaire à coin tronqué 


11.2.4.1 Première étude : comparaison directe 

Nous évaluons à IxS donné en fonction de la distance r d'éloignement, les champs 
magnétiques respectifs de la boucle élémentaire et d'une boucle carrée. Le point 
d'observation est dans le plan des boucles (même circuit imprimé). 

La figure II.17 donne l'amplitude maximale Hy du champ (évaluation par code de 
couleurs) pour les deux boucles, à z = 0. 

Les paramètres pour le tracé sont : f = 100 kHz, Se = 25 ст, I = 10 A. 


Y (m) 
Y (m) 


X (m) 


Ы 060 


Figure П.17: Amplitude maximale du champ magnétique pour les deux boucles 
BCD 061-062 : Influence de l’éloignement (figures П. 17а et II.17b) 


On observe que pour r supérieur à deux fois le côté du carré, les amplitudes sont 
pratiquement équivalentes. On notera également que le champ est maximum au 
centre de la boucle élémentaire, ce qui est normal puisqu'elle est supposée 
« infinitésimale », alors que le champ le plus élevé (d'amplitude plus faible) se situe 
au voisinage du fil pour la boucle carrée. 


П.2.4.2 Deuxième étude ` variation de la taille de la boucle carrée 
Nous proposons maintenant d'examiner l'influence de la taille de la boucle carrée 
sur l'amplitude du champ magnétique en des points situés dans le plan de la boucle. 
Pour des raisons de facilité de représentation, les points suivent une trajectoire 
rectiligne partant du centre de la boucle et perpendiculaire à un cóté (figure II.18). 
Nous proposerons ensuite une variante de l'étude dans un plan paralléle à la 
boucle à quelques millimétres au-dessus de celle-ci. 
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0.03 


Direction 
d'observation 


-0.02 


-0.03 


n 


-0.03 -0.02 -0.01 0 0.01 0.02 0.03 
X 


Figure 11.18 : Représentation des différentes boucles carrées étudiées 


La figure П.19 représente l’évolution de l'amplitude maximale (Hy) du champ 
magnétique en fonction de l'éloignement r, pour les deux types de boucles avec, 
pour la boucle carrée, un paramétrage par la longueur a du cóté. 

L'observation est réalisée en 2D dans le plan des boucles (XY). 


DRE epos eee Se Heec x dd nin ie oe | 
- а=0.005_. : —— Boucle carré | 

Огеш» DAP tegere eS te Boucle élémentaire |` 
gor |- a= 001. DET IRE ec nes ent TETEE ЕБЕ Р 

2 60 ` зорам е 2 * Paramétres : 

a À | ў | i I=10 A m 
ENS "AES a 771 Yp=0 Zp=0 +. 
"BS N JÀ ИИ) 

0 Марр anc 24 e us cad Eat, oe < dr, Leit, ipo оа äisen, A aut 
Юр E LE = Ze Rar | 
$ se 

-40 А í А А 1 
0 0.02 0.04 0.06 0.08 0.1 
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23063 
Figure П.19 : Evolution de Hy pour les différentes boucles 


RCD 064 : Variante en 3D (figure H.19a) 
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Observations essentielles : 

- pour les spires carrées, Hy présente des pics très marqués quand r = a/2 
(situation sur le fil) ; 

- ат = 0, (au centre des spires) les valeurs de Hy sont distribuées en rapport 
avec la pondération (à IXS donné) entre la valeur du courant et celle de la surface, 
ce qui est normal ; 

= les courbes Hy(r) de la boucle carrée se confondent progressivement avec 
Hy(r) de la spire élémentaire pour des valeurs de r élevées. 

A l'issue de ces observations, il semble que l'équivalence sous condition des 
boucles élémentaires (Bg) et carré (Bc) soit envisageable. Pour examiner cette 
possibilité avec un minimum de rigueur nous proposons maintenant d'étudier 
l'erreur relative entre les champs des deux boucles types, soit : 


A% = (Нм )gc — (Au) gg x100 (11.79) 


(H M Jac 
La figure П.20 représente les courbes correspondantes A %(R), paramétrées par 
(a). On observe la convergence asymptotique vers 0 % dès que r est suffisant. 


100, = - > - EE EE 


al... EE 


Erreur relative en % 


0.1 0.2 0.3 04 05 0.6 
r (m) = 965 


Figure 11.20 : Erreur relative sur l'évaluation Hy entre les deux boucles en fonction de r 


La figure IL21 permet d'apprécier l'équivalence en fonction du rapport r/a de 
l'éloignement au cóté de la boucle carrée. Dans ce cas, l'erreur est sensiblement la 
méme quel que soit a. 

En acceptant 10% d'erreur, il suffit de s'éloigner de 2a pour, sans trop 
d'imprécision, assimiler la boucle carrée à une boucle élémentaire et profiter ainsi 
de la simplification importante des expressions mathématiques à implanter. 
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Figure П.21 : Erreur relative sur l'évaluation de l'amplitude maximale entre les deux boucles en f{r/a) 


Comme annoncé au début de ce paragraphe, nous proposons de déplacer le point 
d'observation et sa trajectoire au-dessus de la boucle dans un plan parallèle 
(élévation suivant l'axe OZ) comme indiqué sur la figure II.22. 


Figure I1.22 : Représentation de la direction d'observation et des variables X, et r 


La variable d'éloignement r (coordonnées sphériques) se confondait avec la 
variable X, des coordonnées cartésiennes au point d'observation dans l'étude 
précédente, maintenant, les variations de Hy peuvent être observées en fonction de r 
ou de X, de manière différenciée. 

La figure 1.23 représente les courbes Hu(X,) paramétrées par a pour la boucle 
carrée, l'observation étant réalisée à 1cm au-dessus de ces boucles. 
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Figure 11.23 Evolution de l'amplitude H„ en fonction de X, 
RCD 068 : Variante en 3D (figure П.23а) 
On observe que les pics précédemment évoqués sont maintenant trés amortis, et 
que la convergence vers la courbe Hy(X,) tracée pour la boucle élémentaire se 
produit pour des valeurs de l'éloignement plus réduites. 


La figure 11.24 montre que l'erreur relative est plus faible que précédemment à 
éloignement identique. 
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Figure П.24 : Erreur relative sur l'évaluation de Hy entre les deux boucles en fonction de ғ 
ХС” 070 : Variante en 3D (figure П.24а) 
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Remarque : si au lieu de déplacer le point d’observation parallèlement à la boucle 
au-dessus de celle-ci, on lui fait suivre une trajectoire radiale (8 fixé et r variable), 
l'erreur relative sur Hw, en fonction de l'éloignement, est différente suivant le choix 
initial de 0. On vérifie ainsi que l'équivalence entre les deux types de boucles est 
aussi conditionnée par 0. 


1.2.4.3 Troisième étude ` variation des dimensions à surface fixée 

Toujours dans l'objectif d'examiner l'équivalence entre les deux boucles, nous 
allons maintenant étudier l'influence des dimensions d'une boucle rectangulaire. 

Nous choisissons de fixer le courant et la surface de la boucle dont le rapport 
entre deux cótés adjacents sera néanmoins variable de 0,2 à 5. 

La figure 11.25 illustre ce choix. Le point d'observation se situe dans le plan de la 
boucle et sur un axe passant par le centre perpendiculaire à l'un des cótés. 


| Direction į 
0.01 ALL | : d'observation į 


Y (m) 
T A 


- | 8=0.5 
À -6- 0.33333 
-0.02 Mmm l= 0.25 
| Hi ME | 
-0.02 -0.01 0 0.01 0.02 
X (m) 


Figure 11.25 : Représentation des différentes boucles et de la direction d'observation 


La figure П.26 représente l'évolution de l'amplitude maximale (Hy) en fonction 
de l'éloignement r, avec pour paramétre à. 

- Au point (0,0), les écarts observés sont liés au rapprochement ou non des fils. 
Pour les valeurs extrêmes (0,2 et 5) de 8, l'amplitude Hy est maximale, deux côtés 
étant alors trés rapprochés du centre. 

- Les pics observables sont à nouveau trés marqués. Ils ne dépendent que de la 
valeur de ô (I et S fixés). Ils se produisent toujours quand le point d'observation 
est sur le fil avec une faible contribution globale des autres cótés de la boucle. La 
variation d'amplitude des pics est donc non significative sachant que d'autre part 
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la discrétisation introduite au niveau des calculs numériques ne conduit pas à plus 
de précision si l'on veut conserver un temps de calcul raisonnable. 
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0 0.01 0.02 0.03 0.04 0.05 0.06 
r(m) 
: El 
Figure 11.26 : Évolution de Hy en fonction de l'éloignement r paramétrée par б 


Әср 072 : Variante en 3D (figure П.26а) 


La figure II.27 donne l'erreur relative en fonction de r entre le champ de la boucle 
rectangulaire considérée et celui de la boucle carrée (ô = 1). 
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Figure 11.27 : Erreur relative sur Hy entre les différentes boucles et la boucle carrée en fi r) 
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Les trajectoires observées paramétrées par ô, «encadrent» la valeur 0 


correspondant évidemment à 6 = 1. Ces trajectoires convergent toutes vers zéro dès 
que r est suffisante. 


La figure IL28 représente à nouveau l'erreur mais en fonction du rapport 
« r/longueur / du côté parallèle à l'éloignement ». 


10 


Erreur relative (96) 
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Figure 11.28 : Erreur relative sur Hy entre les différentes boucles et la boucle carrée en f(r/) 


Toutes les courbes convergent franchement vers zéro avec г//. Si l'on tient compte 
du cóté considéré et du sens d'éloignement conjointement, la convergence est trés 
comparable. 

On peut considérer si le rapport longueur/largeur d'un rectangle n'est pas excessif 
ni dans un sens ni dans l'autre que l'équivalence sous conditions est acceptable, 
comme en témoigne à nouveau les isoclines de la figure П.29. 

Avec l'éloignement, toutes les amplitudes s'atténuent et tendent vers zéro, comme 
illustré par les figures complémentaires du CD d'accompagnement. 

En conclusion, nous pouvons par «équivalences » successives, sous certaines 
conditions évoquées auparavant, assimiler une boucle rectangulaire à une boucle 
élémentaire. 
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Figure 11.29 : Courbes d'égale amplitude du champ H pour б=5 
MCD 076-077 : Figures complémentaires (figures II.29a et II.29b) 


1I.2.4.4 Quatrième étude : troncature 


Nous proposons dans cette dernière étude, comme évoqué en début de paragraphe, 
de déterminer les conditions d’équivalence entre une boucle rectangulaire à coin 
tronqué et la boucle élémentaire. Nous ne pouvons bien sûr pas étudier tous les cas 
possibles. Nous retiendrons donc celui de la boucle carrée avec une «enclave » 
égale à 25 96 de la surface (figure II.30). 

De maniére à effectuer une étude d'équivalence avec la boucle élémentaire dans 
les meilleures conditions, il est opportun de centrer cette derniére au centre de 
gravité (Og) du « carré tronqué » et de retenir une direction d'observation respectant 
un minimum de symétrie. 

Aprés modélisation et simulation on obtient l'évolution de l'amplitude Hy en 
fonction de l'éloignement (figure IL31), suivant la direction d'observation (figure 
II.30). On observe que deux pics se produisent correspondant aux valeurs r, et r, de 
r, ils « encadrent » le pic unique du champ de la boucle élémentaire. Ces pics sont 
atténués bien sûr avec l'éloignement, comme illustré sur le CD. 

Les champs de la boucle tronquée et de la boucle élémentaire convergent 
rapidement avec l'éloignement vers la même valeur. 

La figure 11.32 donne l'erreur relative sur le champ de Вт comparé à celui de Bg. 

L'erreur tend à s'annuler à nouveau rapidement avec l'éloignement mais dépend 
fortement de la direction retenue, la plus favorable correspondant au choix 
précédent. 
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Figure II.31 : Evolution de Hy pour Bg et Br suivant la direction d'observation choisie 


WCD 079-080 : Influence de l’éloignement (figures П.31а et IL.31b) 
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Figure 11.32 : Courbes représentant l'erreur relative sur H entre Br et Bg 
ECD 082-083 : Influence de l’éloignement (figures П.32а et II.32b) 


En conclusion, l’équivalence entre les boucles Вт et Bg au regard de l'évolution 
du champ Hy serait encore possible, sous réserve d'ajouter aux conditions 
restrictives des études précédentes une condition liée au respect de l'équilibre de la 
distribution des surfaces de la forme tronquée. Si l'on s'éloigne trop de cette 
condition supplémentaire, comme par exemple dans le cas de la boucle tronquée de 
la figure 1L33 qui présente une partie trés étroite, on peut avoir une erreur 
d'amplitude trés différente en périphérie de la forme, comme en témoigne la figure 
11.34. En effet sur celle-ci, sont représentées les courbes d'égale erreur autour de Br 
témoins de cette derniére remarque. 
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Figure 11.33 : Boucle émettrice très dissymétrique 
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Figure 11.34 : Erreur relative sur H entre la boucle dissymétrique et Bg 
SCD 085-086 : Influence de l'éloignement (figures IL.34a et IL.34b) 


IL.2.5 Limite de validité des équations simplifiées dans l'hypothèse du « champ 
proche » 


Avant d'étudier la validité des équations simplifiées, il semble opportun de 


présenter le paramètre Br et son influence sur l'évolution temporelle du champ 
magnétique. 


П.2.5.1 Présentation du paramètre Вг 


Partant d’une boucle carrée, on peut, dans l’hypothèse d’un éloignement faible, 
apporter un minimum de simplifications aux équations complètes permettant ainsi 
d’obtenir une résolution numérique plus rapide. 

Avec le respect des conditions d’équivalence édictées précédemment, on peut 
raisonner également sur la boucle élémentaire et obtenir en « champ proche » les 
équations beaucoup plus simples présentées auparavant. 

Les résultats d'analyse de validité sont trés comparables, nous exposons 
simplement ici ceux qui concernent la boucle élémentaire. 

Les limites permettant de qualifier de champ « proche » ou « lointain » dépendent 
du paramétre mixte Br (produit nombre d'ondexéloignement). 

Nous jugeons opportun d'étudier préalablement ce produit Br en fonction de la 
fréquence et de l'éloignement (figure II.35). 

La vision 3D (figure II.36) permet d'observer entre autre la « zone de transition » 
entre le «champ proche» et le «champ lointain» et conduire une premiére 
réflexion connexe sur l'influence de la fréquence et/ou de l'éloignement. 
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Figure 11.35 : Evolution de Br en fonction de l'éloignement r et de la fréquence 
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Figure 11.36 : Evolution de Br en fonction de l'éloignement r et de la fréquence en 3D 


Suivant que la valeur de Br est très supérieure, très inférieure ou voisine de 1, on 
sait que l'on est à priori en « champ lointain », en « champ proche » ou en « zone de 
transition ». La représentation en 3D permet dans chaque cas de conduire une 
réflexion sur la pondération entre les valeurs de r et f. 
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II.2.5.2 Influence de Br sur l’évolution temporelle du champ 

Jusqu'à présent nous n'avons retenu que la valeur maximale du champ. Nous 
allons ici observer son amplitude au cours du temps, le courant d'excitation supposé 
sinusoidal de la boucle étant discrétisé en une vingtaine de points. La figure П.37 
donne la représentation de H(r) à l'aide de vecteurs, pour Ө fixé à 45 ° et pour 
différentes valeurs de Br. 

Chaque vecteur correspond à un instant donné. 


Br = 0, 04 Br = 0, 80 Br = 2 
zone de champ proche zone de transition zone de transition i 
unidirectionnalité modification azimutale ^ modification azimutale | 


i zone de transition zone de transition zone de champs couplés: 
: modification azimutale modification azimutale unidirectionnalité 


Paramètres : Ï = 10 A, f = 100 kHz, Se = 25 cm?, Ө = 45 ° 


Figure 11.37 : Evolution de H(t) en fonction de Br 


On peut constater qu’il est trés difficile de trouver une orientation privilégiée en 
zone de transition, illustrant ainsi les difficultés d’analyse des convertisseurs de 
puissance haute fréquence. 

La figure IL.38 représente H(t) à Br fixé et Ө variable de 0 à 90 °. 
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On peut observer que la direction n’est plus unique à chaque instant pour les 
valeurs intermédiaires de 6. Cela traduit le début du passage de la zone de « champ 
proche » à la zone de transition. Le changement minime sur la figure va s’accentuer 
avec l'augmentation de Br, ce qui rend difficile la détermination d'une orientation 
privilégiée en un point de l'espace en zone de transition. 
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Figure 11.38 : Evolution de H(t) à Br fixé en fonction de 0 


11.2.5.3 Calcul du champ avec et sans simplification 

La figure IL.39 représente l’évolution du module de Hy en fonction de Br dans le 
plan de la boucle obtenu par le calcul simplifié en: « champ proche » et le calcul 
exact. On observe une parfaite convergence pour Br < 1, puis une dispersion 
progressive au-delà de 1. 

La figure П.40 donne l'erreur relative sur le champ en un point d'un hémisphère 
au-dessus de la boucle, en fonction de Dr et pour différentes valeurs de 0. 

On observera que l'erreur est rapidement importante dés que Br devient supérieur 
à 1, et que Ө a une influence considérable également (variation de +/— 20 % environ 
pour Dr < 1). Dans ces conditions, il était intéressant d'obtenir la plage d'erreur 
pour Br < 1. la figure П.41 présente les courbes correspondantes, à partir desquelles 
on peut proposer des limites de validité par exemple : 

= sil'ontolére 15 % d'erreur, quel que soit Ө, on prendra Вг < 0,6 ; 

= pour 5 % d'erreur on prendra Вг < 0,3 correspondant àr # Лбл. 
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Figure 11.40 : Erreur relative sur H en fonction de Br pour différentes valeurs de Ө 


EJCD 092 : Variantes en 3D en fonction de 0 (figure П.40а) 


En conclusion, nous retiendrons que l'usage des équations simplifiées est 
envisageable mais conditionné par une valeur de Br « franchement » inférieure à 1. 
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Figure II.41 : Erreur relative sur H pour Br < 1 
EICD 094 : Vue en 3D en fonction de Ө (figure IL41a) 


En optant pour la valeur 0,3 il suffit enfin de vérifier que la direction de H enun 
point quelconque et pour chaque instant est pratiquement unique. 

La figure IL42 tracée pour Br = 0,3 et 0€[0,90 °] illustre cela (un vecteur 
représente H(t) à t donné). 
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Figure 11.42 : Evolution de H(t) à Br = 0,3 pour différentes valeurs de 8 
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11.2.6 Influence d'un plan de masse éventuel 


Comme présenté dans l'introduction générale, les auteurs de ce livre veulent 
proposer une piste de réflexion pour définir une position idéale des cartes sensibles 
dans le boîtier supposé non métallique d'un constructeur. Il est évident que l'intérêt 
est de minimiser ainsi le nombre d'écrans à mettre en place. 

Dans ces conditions édictées par des considérations technico-économiques, la 
présence d'un plan de masse peut étre non souhaitable. C'est pourquoi l'étude qui va 
suivre a été effectuée à priori sans ce plan de masse et que les auteurs se sont 
contentés de montrer rapidement les effets essentiels qui seraient liés à sa présence. 
Ils s'intéressent uniquement au champ magnétique proche. Le plan de masse 
éventuel en cuivre qui serait disposé à l'encontre du champ électrique a une 
influence réduite sur le champ magnétique. Une partie du champ est absorbée et une 
faible partie du champ restant est réfléchie. 

Seule est présentée l'influence de l'écartement « boucle - plan de masse » sur le 
champ résultant aprés absorption et tenant compte d'une réflexion limitée en 
amplitude à 20 % par exemple. Compte tenu de ces hypothèses restrictives, 
l'application de la théorie des images suffira. Voici ainsi l'influence d'un champ 
réfléchi émis par une boucle image de la boucle réelle pour différentes valeurs de 
l'écartement boucle - plan. 

La figure IL43 donne l'amplitude maximale Hwy résultant, dans un plan 
perpendiculaire à la boucle coupant celle-ci en son centre suivant OX et situé à 1 cm 
au-dessus. 
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Figure П.43 : Evolution de Нм en fonction de la hauteur du plan de masse 
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La courbe la plus élevée correspond au champ que l'on obtiendrait sans plan de 
masse. Les autres courbes sont paramétrées par la distance entre la boucle et le plan 
de masse. 

On observe : 

- une réduction du champ quand le plan de masse est trés proche de la boucle 

(face cuivrée du circuit imprimé par exemple) ; 

- deux maxima qui correspondent à l'influence du fil (étude précédente). 

La figure П.44 montre pour Z = 5 cm que l'orientation du champ résultant peut 
étre un peu différente de celle du champ initial. 
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Figure 11.44 : Champ émis avec et sans plan de masse 


IL2.7 Etude simplifiée de l'émission conjuguée de deux boucles 

Les boucles à priori de tailles différentes, sont supposées proches et dans un même 
plan (exemple de la figure II.45). 

Il est trés difficile d'envisager tous les cas de figure suivant les fréquences, les 
amplitudes et les phases des courants dans les deux boucles. La plupart des logiciels 
qui permettent d'observer le champ électrique des cartes numériques supposent que 
les signaux sont identiques dans les différentes pistes excitées. Nous nous 
contenterons ici d'une approche similaire pour le champ magnétique. Les 
paramétres essentiels sont donc géométriques. 

La figure [1.46 donne l'amplitude du champ résultant pour les deux boucles de la 
figure IL45, excitées par le même signal sinusoidal avec Z = 2 cm, Y = 0 et X 
variant de —0,8 à 0,1 m. 

L'orientation dans l'espace du vecteur H résultant est évidemment également un 
peu modifiée. 
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EE ees dad gels Chane Mera 


0.02 рее 


-0.03 а, i П 1 i 


-0.03 -0.02 -0.01 0 0.01 0.02 0.03 0.04 0.05 
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Figure 11.45 : Position des deux boucles émettrices 


dans le plan XY 


iL 


H résultant 
HB2 


H (A/m) 
ч 
2 


Paramètres : 
І= 10 A 
Е = 100 kHz 
Z=2cm 
écart = 4 cm 
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-0.1 -0.05 0 0.05 
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Figure 11.46 : Amplitude des champs émis par chaque boucle et champ résultant 


EICD 099 : Variante pour deux boucles plus éloignées (figure П.46а) 


Les figures П.47 et IL48 donnent les courbes d’égale amplitude du champ 
magnétique dans un plan parallèle aux boucles, respectivement, pour Z = 2 cm et 


0 cm. 
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Chapitre II Présentation simplifiée du rayonnement des cartes de puissance 


[H] (dB) 


i Paramètres : 
i 1= 10А 

i F = 100 kHz 
ï Z=2cm 


X (m C 
(m) Æ 100 
Figure 11.47 : Courbes d'égale amplitude du champ H pour Z = 2 cm 
EACD 101 : Variante pour deux boucles plus éloignées (figure П.47а) 
[H] (dB) 
Paramétres 
I=10A { 
F= 100 kHz | 
Z=0 H 
X (m) | | 
Bio 


Figure 11.48 : Courbes d'égale amplitude du champ H pour Z = 0 
BCD 103 : Variante pour deux boucles plus éloignées (figure II.48a) 
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Chapitre II Présentation simplifiée du rayonnement des cartes de puissance 


En conclusion, l’action conjuguée des deux boucles est un problème délicat à 
traiter dans le cas général, d'autant plus que la proximité des fils modifie fortement 
l'amplitude du champ résultant comme en témoignent les courbes tracées dans le cas 
simple retenu. 

Dans un contexte de CEM, où l’on estime souvent les perturbations extrêmes, on 
peut imaginer que l'on ajoute directement les champs en supposant qu'ils ne sont 
pas déphasés (certains logiciels du commerce procédent de cette facon). 


Conclusion du chapitre 


Dans ce chapitre, nous avons présenté la notion de « boucle équivalente » liée à 
une « topologie » d'un convertisseur de puissance. Nous avons ensuite estimé les 
limites raisonnables fixant les conditions d'équivalence possible entre cette boucle et 
la boucle dite élémentaire. Les paramétres étant liés à l'éloignement, la fréquence, la 
taille de la boucle et les dimensions du conducteur correspondant. 

Au travers de multiples illustrations, nous avons montré comment vérifier que les 
facteurs déterminants sont essentiellement le produit Br et le rapport entre les 
dimensions de la boucle et l'éloignement. 

Les boucles ne sont pas en général de simples rectangles, c'est la raison pour 
laquelle nous avons vu le cas particulier d'un carré « tronqué ». L'influence de cette 
modification géométrique a été examinée en comparant le champ émis à celui de la 
boucle élémentaire. 

Les dimensions du conducteur ont été également prises en compte dans une étude 
comparative s'appuyant sur le calcul direct et l'application de la méthode des 
éléments finis. 

Enfin, nous avons évoqué l'influence possible d'un plan de masse et l'action 
conjuguée de deux boucles coplanaires. 


En résumé, l'ensemble de ces investigations a permis, sous certaines conditions 
explicitées, de dégager la notion de boucle «équivalente » pour la recherche 
simplifiée du champ magnétique proche émis par une carte de puissance. 

Dans le chapitre suivant, où nous allons étudier le signal reçu par une boucle 
victime, nous allons retenir cette notion de «boucle équivalente » de forme 
simplement circulaire ou carrée supposée réalisée avec un conducteur de section 
négligeable. 
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CHAPITRE III 


F.E.M. INDUITE 
DANS UNE BOUCLE SENSIBLE 
ET 
IMMUNITÉ 


La démarche présentée ici est limitée à la recherche d’une stratégie de réduction 
de la tension que peut induire une boucle perturbatrice (Bp) dans une boucle (B,) 
victime voisine, réceptrice du champ H. 

B, est une boucle topologique que l'on assimilera à un simple carré ou à une 
boucle « équivalente » comme explicité au chapitre II. 

La réduction de la tension recueilie par B, sera simplement obtenue ici par 
optimisation de son positionnement par rapport à Bp. Avec cette démarche, on 
évitera l'utilisation d'écrans. 

L'étude est conduite sur une boucle Bp unique correspondant à un « ilot » tel que 
présenté au chapitre II. B, est une partie de la carte de commande/contróle d'une 
structure de puissance incluant Bp, ou appartient à une carte additionnelle de 
traitement de signal d'un capteur extérieur (température par exemple) associé au 
convertisseur. 

La boucle réceptrice B, peut ainsi étre constituée par les connexions entre la sortie 
d'impédance Zs et l'entrée d'un autre circuit (d'impédance Ze). Sa géométrie 
dépend de la trajectoire des pistes (liaison entre noeuds et liaisons au potentiel de 
référence, 0 V ou 4 Vcc). 

L'impédance de cette boucle potentiellement victime est donc liée à la technologie 
des circuits et aux longueurs des pistes. 

La recherche du courant dans la boucle victime est difficile car liée à de nombreux 
facteurs (technologiques en particulier). D'autre part, nous savons que les tests 
habituels d'intégrité de signal procèdent par «excitation » en tension d'une piste. 
Cela nous a conduit à retenir essentiellement la f.e.m. induite comme grandeur 
perturbatrice, pouvant introduire par exemple dans le cas d'un amplificateur 
opérationnel un bruit préjudiciable au signal utile. 

Comme précédemment, nous avons eu le souci de simplifier les écritures à 
implanter, et avons en conséquence conduit une étude comparative tenant compte ou 
non des simplifications introduites. 


Le chapitre se décompose ainsi : 

- étude théorique pour la recherche de la f.e.m. induite e(t) incluant plusieurs 
démarches simplificatrices comparées ; 

- cartographie de e(t), les boucles B, et B, étant dans des plans 
perpendiculaires ou paralléles correspondant à des assemblages usuels industriels ; 
- recherche de la position optimale de B, tenant compte des contraintes 
techniques de fixation de la carte sensible. 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


III.1 Recherche et estimation de la f.e.m. induite 


Dans le souci de simplification souhaitée, plusieurs démarches sont entreprises 
pour la détermination de e(t). Pour les présenter de manière cohérente, nous 
proposons un organigramme incluant des blocs de calcul repérés par les chiffres 1 à 
5. Les calculs seront explicités et les résultats comparés. 


III.1.1 Organigramme de la recherche de e(t) 

La représentation est donnée figures III.1 et Ш.2. 
Description générale 

- X0) est le courant d'entrée dans la boucle perturbatrice Bp, il est choisi de 

forme trapézoidale tel que défini au chapitre I avec ou sans oscillation due par 

exemple à un CALC associé au composant en commutation ; 

-  i(f)est ce même courant dans le domaine fréquentiel (spectre) ; 

- e(t est l'onde temporelle de flux ; 

- e(t) est l'onde temporelle de la f.e.m. induite dépendante de X, Y, Z ou 

X,Y,Za; 

- = e (D), ef), es(f), e4(f) et es(f) sont les spectres de la f.e.m. induite déterminés 

pour les blocs de calculs correspondants. 

Concernant le paramétrage [X,Y,Z,a] il convient de préciser que la boucle 
réceptrice B, est centrée au point P de coordonnées Xp, Үр, Zp (Bp étant centrée en 
[0,0,0]) et en rotation possible (repérée par un angle œ) autour d'un diamètre 
parallèle à OX (figure III.3). 

-  uS,ii.H représente le produit scalaire du vecteur normal à la surface de В, par 


le champ magnétique H multiplié par uS, (produit perméabilité x surface de B,) 
de manière à déterminer le flux $(f) ; 

-  y$,H représente le produit de la surface de B, par la valeur maximale (notée 
H) du champ au point Xp, Yp, Zp en supposant le champ normal à B,. Ce qui 
correspond également au flux maximal que peut recevoir B, (pire des cas au sens 
de la CEM) ; 

- la notation de bloc « F.F.T. » signifie que l'on applique la transformée de 
Fourier rapide au signal d'entrée (calcul numérique) «LF.F.T.» est Ја 
transformation inverse ; 

- la notation de bloc «dérivation par rang» signifie que l'on applique 
l'opération de dérivation (par rapport au temps) à chaque onde sinusoidale de 
fréquence harmonique constituant le spectre du flux ; 

- la notation « dérivation » signifie que l'on effectue directement la dérivation 
du flux ф(0 ; 

- la question «champ uniforme?» signifie que l'on s'interroge sur la 
possibilité de considérer qu'en tout point de B, , le champ H а la méme amplitude 
et la même orientation. Les figures III.4 et IIL5 illustrent cette démarche ; 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


CALCUL 
Recherche de di 


(D 


Pire des cas? 


Dérivation 


e(t) ou ext) 
paramétrage par 
X, Y, Z, еа 


(1) : Bloc de calcul général 


(2) : Bloc de calcul avec l'hypothése du champ uniforme 


Figure Ш.1 : Organigramme (1/2) 


81 


// 


Chapitre III F.e.m. induite dans une boucle sensible et immunité 


Simplification 


RAS Non 
des équations ? 


Oui 


Boucle 


élémentaire 
(4) 


Estimation par la 
méthode des é( f)=uS,H 


enveloppes ? (3) 


Oui 


Calcul des points 
particuliers (5) 


es(f) : 


enveloppe de 


ea(f) 


Dérivation 


ext) 
paramétrage par 
X,Y et Z 


ext) 
paramétrage par 
X,YetZ 


(3) : Bloc de calcul avec les hypothèses du champ uniforme 
et du flux recueilli maximum 

(4) : Bloc de calcul avec le maximum de simplification 

(5) : Bloc de calcul basé sur la notion d'enveloppe 


Figure Ш.2 : Organigramme (2/2) 
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Chapitre ITI F.e.m. induite dans une boucle sensible et immunité 


-  laquestion « simplification des équations ? » signifie que l'on s'interroge sur 
la possibilité d'appliquer les relations de la boucle élémentaire équivalente en 
champ proche à la place des relations générales de la boucle carrée retenue 
initialement ; 

~ la question «estimation par la méthode des enveloppes ? » concerne la 
possibilité de retenir essentiellement les points de l'enveloppe des maxima des 
harmoniques du spectre de i(t), correspondant aux changements de pentes. 


Axe de rotation 


S + 
Vue en coupe (YZ) pour 9 z 0 
Figure III.3 : Position et repérage de Br et Bp 
0.04 
ml _: | | = fer 
0.02 | - Ae 
єз 0.01 X 
0 S - 2 ї 
d К 
0.01 З 21 
-0.02 > AE ll 
X (m 2 0.08 
(m) 0.0 Eon 0 0.02 0.04 0.06 
m Æ 104 


Figure III.4 : Champ quasi uniforme sur la boucle B, 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


Figure Ш.5 : Champ non uniforme sur B, 


III.1.2 Blocs de calcul 

Le calcul s'effectue dans le domaine fréquentiel en considérant que les 
harmoniques de champ ont pour amplitude maximale leur valeur à t = O. Cela limite 
la validité du calcul à un éloignement raisonnable. Par exemple, pour Br < 0,3, 
l'erreur relative maximale sur l'amplitude du courant est inférieure à 5 96. 


III.1.2.1 Bloc de calcul général 

Le calcul du flux $(f) peut être effectué sans hypothèse d'uniformité du champ, 
avec les équations exactes pour le champ émis par la boucle carrée dont toutefois on 
négligera la section du conducteur et les effets d'interaction avec une boucle 
réceptrice équivalente circulaire. 


Hf = [JEG ds, = Јоха) ds, = jac dl (Ш.1) 
5 S, І 


r r 


Nous pouvons réaliser le calcul de l'intégrale curviligne par la méthode des 
trapèzes. Il suffit pour cela de discrétiser en N points le contour de la boucle 
réceptrice comme indiqué sur la figure III.6. 

Pour chacun de ces points, nous avons déterminé le vecteur potentiel At f) à 


partir des équations du chapitre précédent et réalisé le calcul du produit scalaire. 
À;(f) di 
Pour -90 < œ < 90: dl; =—sin(y; )x'+ cos(y; )y' (II.2) 
avec: у'= соѕ(о) у tsin(a)z et x'- х 
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Chapitre ITI F.e.m. induite dans une boucle sensible et immunité 


Donc : À;(f)- dli = A Cf)sin(y;) + A,Cf)cos(y;)cos(a) ШЕ 


Figure III.6 : Représentation de la discrétisation de la boucle réceptrice dans son plan 


Le contour de la spire réceptrice étant circulaire, l’expression du flux est donc la 
suivante : 


^ . A i A Á. dl, AC dl 
600 = Län A =Y. rx(yi - BEER IT ша 
1, i=2 


III.1.2.2 Bloc de calcul avec l'hypothèse du champ uniforme 

Une première simplification peut être introduite avec l’hypothèse de l'uniformité 
du champ recueilli par B, mais l'utilisation du produit scalaire Я.Н permet toujours 
de tenir compte de sa position (Xp, Yp, Zp) et de son orientation (œ). 

Bien entendu l'hypothése de l'uniformité du champ a une limite de validité liée au 
rapprochement des boucles. 

Une remarque à ce sujet sera présentée au paragraphe (III. 1.3). 

Calcul du flux à partir du champ au centre de la boucle : 
ф= [8.ás, = ugS, H, i 

S 

ñ est la normale à la surface de la boucle réceptrice et Н le champ magnétique au 


(Ш.5) 


r 


centre de celle-ci. 
L'angle o entre H с et le plan de la boucle est déterminé en exprimant H. dans la 


base (x',y',z'). M est alors la matrice de passage des coordonnées de la base 
(X, y, z) à la base (x', y',z"). 
1 0 0 
M =|0 cos(œa) sin(a) 
О -sin(œa) cos(æ) 
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Chapitre HI F.e.m. induite dans une boucle sensible et immunité 


donc : 


H 
Q'— arctan c Z 01.6) 


le calcul du flux est donc le suivant : 


T o. di nl: 
9 = uo, H. cos =a ) = yos, Ho sin(a’) (III. 7) 


III.1.2.3 Bloc de calcul avec les hypothèses du champ uniforme et du flux recueilli 
maximum 

Cette hypothése constitue une étude dans le « pire des cas », correspondant donc 
au maximum possible de f.e.m. induite obtenu lorsque le champ est perpendiculaire 
a B,. Le flux est donc calculé dans le domaine fréquentiel avec les harmoniques 
d'amplitude maximum. 

Dans ces conditions le flux s'exprime simplement par : 

9 = ug S, H. (111.8) 


III.1.2.4 Bloc de calcul avec le maximum de simplification 
Nous mettons à profit ici la notion de boucle équivalente présentée au chapitre II 
en assimilant la boucle carrée à une boucle élémentaire centrée en [0,0,0] dans le 
plan X,Y,Z. Par ailleurs, l'hypothèse du champ proche étant retenue, les équations 
sont très simplifiées dés le départ et autorisent l'obtention directe de la f em. : e4(f). 
Le flux est toujours donné par ф = uoS,H o- 


H, est déterminé à partir des équations simplifiées de la boucle élémentaire en 
champ proche donc : 


IS, [3cos(20)+5 


H,=—& (Ш.9) 
^^ олу? 8 
La tension induite s’exprime en régime sinusoïdal de la façon suivante : 
dy d 
e(f)= = ~ 7, WoHes,)= — jwugH,S, (Ш.10) 
d’où 


3cos(28 )+5 


А IS 
e4(f) = - jwuoS, —5 | (Ш.11) 
2лғ 8 


III.1.2.5 Bloc de calcul basé sur la notion d'enveloppe 


Il s'agit ici essentiellement de retenir les points de changement de pente des 
segments constituant l'enveloppe des maxima du spectre e4(f) déterminé par le bloc 
4. Cette démarche permet d'obtenir es(t) limitée à l'enveloppe de e,(f) résultant de 
l'application sur Bp du signal i(t) toujours supposé trapézoïdal mais sans oscillation 
de commutation superposée. 

Les enveloppes sont présentées dans le paragraphe III 1.3. 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


III.1.3 Application des hypothèses et des méthodes 


III.1.3.1 Calcul général 

Courant i(t) injecté : le courant est de type trapézoidal symétrique (1.3.3.2) avec 
superposition d'une oscillation parasite (équation I.13), dont les paramètres sont les 
suivants : 

- courant: Jus 10 A, f = 100 kHz, к= = 0,1 us 

tu = 2 Us, fp 10 MHz, D,= 10 A; 
- paramètres de B, : centrage en (0,0,0), côté de 5cm, surface de 25 cm? ; 
- Je point d'observation est repéré раг: Ө = 60 °, @ = 90 °, r = 0,05 m. 


© 


Courant I (A) 
an oo 


4 

2 

0 
0 l 2 
temps (s) x 10-5 


Hin 
Figure Ш.7 : Courant parcourant B, (f; = 10 MHz) 
ECD 107 : Courant parcourant B, avec f; = 7 MHz (figure III. 7a) 
Le spectre du courant : i(f) est représenté figure Ш.8. 
On pourra observer que ce spectre est identique à un coefficient peu variable près 
à H(f) et &(. 
Le spectre de la f.e.m. induite, ef, est présenté figure II.9. 
Les paramétres de la boucle B, ainsi que le point d'observation sont inchangés. 
Les paramétres de B, sont : 
- cercle de 2 cm de diamètre centré en P (0 = 60 °, ф = 90 °, r = 0,05 т); 
- @=91° (choisi de manière à obtenir une f.e.m. quasimaximale). 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


Amplitude (dB) 


10° 10° 
Fréquence (Hz) 
NEIS = 108 
Figure Ш.8 : Spectre i(f) du courant parcourant la boucle Bp avec f, = 10 MHz 
ECD 109 : Spectre du courant percourant Bp avec f, = 7 MHz (figure Ш.8а) 
3 E CE TEST 5 md 
AE Sete a sets 
ASP eur „Чү ASE 
са 
3 
3-50} 
< 
-55 
-60 
-65 
10° 
Fréquence (Hz) 
= ] 10 


Figure Ш.9 : Spectre ef) de la f.e.m. induite dans B, avec f, = 10 MHz 
MCD 111 : Spectre ef de la f.e.m. induite dans В, avec f; = 7 MHz (figure Ш.9а) 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


A partir du spectre e,(f) il est possible de reconstituer l'évolution temporelle de la 
f.e.m. induite e;(t). Le tracé correspondant est donné figure Ш.10. 


0.6 


0.4 


c 
FA 


© 


, 
c 
сә 


Tension induite (V) 


-0.4 


-0.8 L J, 
1 2 
temps (s) x 10? 


L3 
5112 
Figure Ш.10 : Evolution temporelle de la Cem. induite еқ) 

Әср 113 : Evolution temporelle de la f.e.m. induite dans B, avec f; = 7 MHz (figure Ш.10а) 


& Remarque concernant les limites du calcul 

Un essai avec un courant i(t) trapézoïdal sans oscillation superposée est réalisé. 

La figure IIL11 présente sa décomposition spectrale et la figure Ш.12 la 
décomposition spectrale du champ magnétique qui lui est associé au point 
d'observation choisi, avec le paramétrage suivant : 

- courant: Iy = 10A, f 2 1 MHz, t, = t; = 0,1 us, ty = 0,2 Us ; 

- В}: сепітаре : (0,0,0), côté 5 cm, surface : 25 cm? ; 

- point d'observation: Ө = 45 °, ф= 0°, г= 1 m. 

On observera des variations d'amplitude rapides liées à des conditions 

paradoxales d'étude. En effet, l'excursion poussée en fréquence nous entraîne en 

zone de transition où le maximum du champ n'est plus atteint à t = 0 comme 
retenu en hypothèse de calcul. 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


20 pa —— ————— 


DO LPL: - - 


-40 HA - -|; 


Amplitude (dB) 


-60 H' - DISK WE 8 


-80 ' ti- Au h 


-100 4 
10° 


Eu 


Fréquence (Hz) 
Figure III.11 : Décomposition spectrale du courant 


#5 —— ——— — — 
1a ` a a ык КА 


Amplitude (dB) 


-130 H >| - H, 


-140 Ш 


10* 


Su 
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Figure Ш.12 : Spectre du champ magnétique H correspondant 


III.1.3.2 Hypothèse du champ uniforme 
Le champ supposé uniforme est calculé au centre de B,. Il a pour paramétres : 
- courant: Iy = 10 A. f = 100 KHz, t = t, — 0,1 Us, ty = 2 us, fp = 10 MHz, 
D,= 10 A; 
- paramètres de B, : centrage (0,0,0), côté 5 cm, surface : 25 cm? ; 
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Chapitre Ш F.e.m. induite dans une boucle sensible et immunité 


- paramètres de B,: cercle de 2 cm de diamètre, centre P (8 = 60 °, @ = 90°, 
r = 0,05 т, œ = 91 °, choisi de manière à obtenir une f.e.m. quasimaximale). 
Le spectre de la f.e.m. induite e,(f) est présenté figure III.13. 


4S5p e 3 2] 
ЕР ie a 
S a 
2 
= D 
OPE Roses spams rene ЫП | 
= 
E 
2 

ont. | IM 

-65 y | | H | | | | 

10° 10° 
Fréquence (Hz) LJ 
Ы 


Figure Ш.13 : Spectre de la f.e.m. induite calculé avec l'hypothèse du champ uniforme 


La f.e.m. induite e;(t) est représentée figure III.14. 


Tension induite (V) 


„—4 


0 1 2 
Temps (s 5 
ps (s) x10 Hu: 


Figure IL 14 : f.e.m. induite calculée avec l'hypothèse du champ uniforme 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


III.1.3.3 Hypothèse du flux recueilli maximum 
Le flux est maximum car le champ en P supposé uniforme est normal à B.. 
Le spectre de la f.e.m. induite ef), (figure Ш.15), a pour paramétrage : 
- 1Ім= 10А, #= 100 kHz, t=t,= 0,1 us, ty = 2 us, fp = 10 MHz Dp= 10A ; 
-  B,:centrage еп (0,0,0), côté 5 cm, surface ` 25 cm? ; 
-  B.:cercle de 2 cm de diamètre, centre P (0 = 60 °, ф = 90 °, r = 0,05 m). 


Tension induite (dB) 


Fréquence (Hz) 


118 
Figure IlI.15 : Spectre de la f.e.m. induite avec l'hypothèse du flux maximum 


f.e.m. induite e3(t) correspondante (figure III.16) : 
0.3 
0.2 


0.1 


Tension induite (V) 


-0.2 


-0.3 


Temps (s) x 10° а, 
Figure Ш.16 : f.e.m. induite avec l'hypothèse du flux maximum 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


III.1.3.4 Hypothèse de la boucle équivalente 


On considère la boucle équivalente (paramètres IxS,) et les équations simplifiées 
correspondant au maximum, dans l'hypothése du champ proche. Le spectre de la 
f.e.m. induite e4(f) est donné figure Ш.17 (même paramétrage que précédemment). 


rrit T T T Y T T ttr T GATI N XX: 


Tension induite (dB) 


Fréquence (Hz) 


= 
120 
Figure III.17 : Spectre de la Lem. induite 
La f.e.m. induite e4(t) correspondante est représentée figure Ш.18. 
= 
= 
ke 
5 
= 
E 
5 
E 
0 1 2 
-5 
Temps (s) x 10 B, 


Figure Ш.18 : F.e.m. induite e(t) 
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Chapitre III F.e.m. induite dans une boucle sensible et immunité 


III.1.3.5 Méthode de l'enveloppe 
On notera les points P1, P2, P3 correspondant aux changements de pentes. 
Nous allons chercher l'enveloppe de i(f) avec les mêmes paramètres 
qu'auparavant, et œ = 91 ? (choisi de manière à obtenir une f.e.m. quasimaximale) 
Le courant est représenté dans les deux domaines figures III.19 et 20. 
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Figure 111.20 : Décomposition spectrale du courant et points caractéristiques 
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L'enveloppe de e(f) est donnée par la figure IIL21. 
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Figure Ш.21 : Enveloppe de la f.e.m. induite 


Remarque concernant les limites de la méthode de l’enveloppe 

Les figures Ш.22 et IIL23 représentent le spectre de la f.e.m. induite pour un 
courant i(t) de forme rectangulaire (signal numérique) et de fréquence élevée 
(21 MHz) à une distance de 10 cm, puis 3 m. 

On constatera qu'il est difficile dans le deuxiéme cas de situer l'enveloppe. Cela 
est bien sür lié à la non-validité des équations, puisqu'à cette fréquence et à cette 
distance (Dr élevé) l'hypothése du champ proche est erronée ! 

Pour que la démarche reste valide il faut toujours que les conditions du champ 
proche soient réunies. 


95 


Chapitre III F.e.m. induite dans une boucle sensible et immunité 


-100 x T "a T 
SE 
1 D D SG Y Xy = H 
= Rx dedu | | 
3-120 nu оти: REIS S EU rl E e Se 
E , Qoa saa í te 
> 
Б f bx А 
56-130 = FOF ët LI р t 
e nr 
* Harmoniques IE E е dE ei 
~ Enveloppe RSC : Ey Reni rni as 
-150 == с = E e x аа a al. A la 
10 10° 10° 
Fréquence (Hz 
q (Hz) Ы 55 
Figure 11.22 : Spectre de la Cem. induite pour r = 10 cm 
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Figure !IL23 : Spectre de la f.e.m. induite pour r = 3 m 
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IIL.1.4 Comparaison des résultats ob tenus 

Les essais avec les blocs 1 à 4 ont été effectués avec superposition d'une 
oscillation parasite ; dans un but comparatif nous proposons d'examiner e(f) à i(t) 
donné. Nous observons donc (figure III.24) conjointement e,(f), e2(f), ex(f) et e(f). 
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Figure DI 24 : Décomposition spectrale de la f.e.m. pour les différents blocs de calcul 


De manière à effectuer la comparaison avec le calcul n°5, nous avons tracé les 
courbes précédentes avec un courant i(t) sans oscillation haute fréquence. La 
comparaison s'effectue en examinant directement les courbes de la figure Ш.25. 
Interprétation des résultats 

Dans la mesure où l'angle o (calculs 1 et 2) est choisi de manière à obtenir un 
niveau induit maximum, il est normal, compte tenu des dimensions des boucles et de 
l'éloignement choisis, d'obtenir des niveaux comparables pour e;(f) et e2(f). 

Si le flux est supposé maximal, il faut s'attendre à un niveau supérieur, pour la 
f.e.m. induite e;(f). 

La f.e.m. e4(f) est trés semblable à e;(f) dans la mesure où les conditions de champ 
proche et d'équivalence de la boucle sont alors réunies. 

On vérifie enfin que e;(f) présente des maxima autorisant le repérage aisé des 
points particuliers permettant d'obtenir sans difficulté l'enveloppe e;(f). 

L'équivalence des résultats est témoin de la validité des calculs et hypothèses. 
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Figure Ш.25 : Décomposition spectrale de la f.e.m. pour les différents blocs de calcul 
En conclusion, il est possible d'évaluer sans trop d'imprécision la f.e.m. induite 
dans une boucle victime B, avec la notion de boucle équivalente pour la boucle 


perturbatrice Bp et les équations du champ proche, sous réserve de respecter les 
conditions rappelées ci-après : 
- dimensions de B, relativement petites devant celles de Bp et devant 
l'éloignement ; 
- forme de B, proche du carré comme explicité au chapitre II ; 
- produit Br (proportionnel à fxr) inférieur à 0,3 ; 
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- position de B, proche de B, mais suffisamment écartée pour éviter les 

champs à proximité des conducteurs réalisant B,. 

Cette évaluation donne bien entendu les valeurs extrémes de la f.e.m. induite (pire 
des cas). 


III.2 Cartographie en situations typiques 


L'étude précédente avait pour but de montrer les possibilités de simplifications 
diverses liées à des conditions précises. 

Dans ce nouveau paragraphe nous retiendrons essentiellement le calcul général 
(bloc 1), car notre objectif est maintenant d'évaluer la f.e.m. induite dans la boucle 
B, dont l'orientation est fixée à priori perpendiculairement ou parallélement à Bp 
rendant impossible le choix particulier précédent de œ qui correspondait à la 
réception maximale du flux. 

De maniére à présenter le plus simplement possible le probléme, nous raisonnons 
sur une seule fréquence pour i(t) laquelle pourrait étre par exemple celle de 
l'oscillation parasite puisque, à priori, la plus critique au regard de la perturbation 
engendrée. 

Le calcul, en particulier la dérivation, est alors trés simple. 


III.2.1 Cas des boucles parallèles 

Dans cette hypothèse, les boucles Bp et B, appartiennent à la même carte de circuit 
imprimé ou à deux cartes distinctes, disposées pour des raisons techniques 
d'assemblage sur des glissières parallèles. 

L'angle œ vaut donc 0° et les seuls paramètres à considérer sont Xp, Yp, Zp 
(coordonnées du centre P de B, dans l'espace). 

La figure III.26 montre une disposition respectant l'hypothèse annoncée. 
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Figure Ш.26 : Disposition possible de B, par rapport à Bp 
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Le champ est représenté uniquement au centre de la boucle sur la figure 
précédente, de maniére à repérer aisément ses changements d'orientation édictés par 
la translation de B, dans le plan parallèle à Bp. 

On voit de suite l'influence significative de la position sur l'amplitude. 

Nous proposons maintenant d'examiner les équipotentielles de la f.e.m induite e(t) 
obtenues par le calcul général (1) pour différentes positions de B, par rapport à Bp. 

Chaque point des équipotentielles est repéré par le point P, centre de B,, de 
coordonnées: Xp, Yp, Zp. Les dimensions des boucles sont fixées comme 
précédemment (5 cm de cóté pour le « carré » Bp et 2 cm de diamétre pour B,). 


Ш.2.1.1 Première position : B, est dans un plan parallèle au-dessus de Вр 
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Figure Ш.27 : Equipotentielles de la Cem induite quand B, est dans le plan de B, (Z, = 0) 


ICD 131-134 : Représentations pour Z, variable de 0,5 à 5 cm (figures Ш.27а à IIL27d) 
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III.2.1.2 Deuxième position ` vues suivant une coupe parallèle au plan YZ 
La figure III.28 présente les équipotentielles correspondantes de la f.e.m. induite. 


Xp =0 Tension induite (V) 


Yp (m CH 
SS 135 
Figure Ш.28 : Equipotentielles de Cem. induite dans le plan YZ pour Xp = 0 


ВСР” 136-139 : Représentations pour X, variable jusqu'à 5 cm (figures Ш.28а à III.28d) 


A noter que l'une des figures du CD correspond au cas particulier intéressant oü 
un fil de coté est situé dans le plan observé (figure III.28b). 

Ces figures bien que non exhaustives permettent de tirer les conclusions 
suivantes : 

- prés des conducteurs, la f.e.m. est très élevée, il faudra donc éviter ces 

positions ; 

- Гахе de la boucle est bien sûr à éviter également ; 

- des zones que nous qualifierons de « tranquilles » notées Zr correspondant à 

des équipotentielles de niveau réduit sont parfaitement repérables et utilisables 

pour positionner B, « au mieux ». 


III.2.1.3 Equipotentielles de la f.e.m. induite dans l'espace autour de Bp 

Le paramétre est la valeur efficace E de la f.e.m. induite qui peut étre édictée par 
le seuil de sensibilité d'un circuit sensible (boucle B,). 

Les figures Ш.29 et Ш.30 de représentation en 3D confirment les remarques 
précédentes. En effet, les zones « tranquilles » se situent à l'extérieur des volumes 
maillés avec la possibilité de se rapprocher suffisamment de la périphérie de Bp avec 
un niveau de f.e.m. restreint. 

Les volumes sont délimités par les équipotentielles de f.e.m. induite. 
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Figure Ш.29 : Equipotentielles de la f.e.m. induite pour E = 15 mV 
BCD 141-144 : Représentations pour plusieurs valeurs de E dégressives (figures Ш.29а à IIL29d) 
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Figure 111.30 : Equipotentielles de la f.e.m. induite dans l'espace pour quelques valeurs de E 


MCD 146-149 : Variantes de représentations dans l'espace (figures III.30a à Ш.304) 


102 


Chapitre III F.e.m. induite dans une boucle sensible et immunité 


III.2.2 Cas des boucles perpendiculaires 

C'est le cas par exemple des cartes en façade d'un appareil perturbées par une 
carte de puissance sur glissière perpendiculaire à la façade. 

La figure Ш.31 montre une disposition correspondante, avec deux positions 
possibles de B, et indication du champ en son centre. On observe le fort changement 
d'orientation de celui-ci. 
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Figure Ш.31 : Disposition possible de В, perpendiculaire à В, 


Nous allons à nouveau examiner les équipotentielles de la f.e.m. induite dans B, 
avec les mêmes hypothèses que précédemment sauf pour l'angle o qui vaut 
maintenant 90 ? (au lieu de 0 ?). 

Deux positions sont prises pour exemple : 

a) un diamètre de B, est dans le plan de Bp (figure III.32). 

Les équipotentielles obtenues sont données figure 1.33. 

On observe deux zones « tranquilles » de chaque côté de la boucle (la position 
intérieure n'étant pas envisageable). 

b) cette fois, le centre de la boucle a pour coordonnées : 

Xp =0, Yp et Zp variables. 

La figure Ш.34 est le tracé des équipotentielles de f.e.m. induite dans le plan YZ, 
Br étant comme annoncé, perpendiculaire à Bp. D'autres tracés disponibles sur le 
CD, permettent de déduire l'influence d'une augmentation progressive de Xp. 

On pourra observer que sur la figure IIL34 sont indiquées 242 zones 
«tranquilles » à l'extérieur de Bp correspondant aux deux cas suivants : Zp # 0 avec 
Yp variable mais supérieur aux dimensions de Bp et Yp # 0 avec Zp variable mais 
supérieur à un minimum (dimension de B,). 
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Plan parallèle 
au plan de B, 


Figure Ш.32 : Le diamètre de B, est dans le plan de Bp 
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Figure Ш.33 : Equipotentielles de la f.e.m. induite dans le plan XY, B, perpendiculaire à Bp(Zp= Іст) 


ВСр 152-153 : Réduction de niveau avec l'accroissement de Zp (figures III.33a HI.33b) 
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Figure III.34 : Equipotentielles de la Cem. induite dans le plan YZ pour B, perpendiculaire à B,(Xp=0) 


ECD 155-156 : Equipotentielles de la f.e.m. dans le plan YZ : Influence de X, (figures Ш.34а et IIL34b) 
ROD 157-160 : Equipotentielles de la f.e.m. dans le plan XZ : Influence de Y, (figures III.34c à Ш.340) 


L'étude de ces deux positions de B, dans un plan perpendiculaire à B, montrent 
que des zones « tranquilles » privilégiées de positionnement de B, sont ici à nouveau 


présentes. 
La figure III.35 


de représentation en 3D confirme cette remarque. En effet, les 


zones « tranquilles » se situent à l'extérieur des volumes maillés avec la possibilité 
de se rapprocher suffisamment de Bp avec un niveau de f.e.m. induite restreint. 
Bien entendu, cette étude n'est pas exhaustive. 
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Figure Ш.35 : Equipotentielles en 3D de la f.e.m. induite pour E = 3 mV 


BCD 161-166 : Equipotentielles en 3D de la f.e.m. induite pour les valeurs dégressives de E suivantes : 
40-30-20-15-10 et 5 mV (figures III.35a à III.35f) 


MCD 167-169 : Equipotentielles de la f.e.m. dans l'espace pour quelques valeurs de E (figures IIL35g à 
IIL35i) 

En conclusion de ce paragraphe, il apparait que dans les cas typiques 
d'assemblage industriel des cartes, il soit possible par simple modification de 
sérigraphie autorisant la «translation » de B, sur le circuit imprimé de la carte 
sensible, de réduire la f.e.m. induite par le champ perturbateur de Bp. 

Pour réduire davantage la f.e.m. induite, il faudrait introduire un degré de liberté 
supplémentaire en se dégageant des astreintes liées aux dispositions techniques 
habituelles. П apparaît néanmoins difficile de faire accepter à des assembleurs des 
positions totalement quelconques conduisant à la réalisation de fixations de formes 
non « d’équerre ». Dans ces conditions, nous n'envisageons comme dans l'étude 
précédente qu'une simple rotation de B, suivant un diamétre paralléle à (OX), 
direction qui peut correspondre à une arête du boîtier non métallique du dispositif 
global. 

Le paramétrage s'effectuera donc en Xo, Yp, Zp et avec œ variable de 0 ° à 90 °. 

En un point donné le positionnement pourra donc être optimisé par œ Le 
paragraphe suivant concerne cette démarche. 
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IIL3 Optimisation du positionn ement de la boucle sensible 


Nous allons dans un premier temps présenter une recherche des valeurs de œ 
permettant de réduire la f.e.m. induite, à Xp, Yp, Zp donnés, cela avec les relations 
générales. Dans un deuxiéme temps nous examinerons les possibilités de considérer 
les démarches de calcul simplifié présentées précédemment. 

Notre démarche se veut pragmatique, c'est ainsi que nous partons d'un cas 
concret. Pour une carte incluant Bp donnée, il existe un boitier réel dont les 
dimensions sont en rapport avec celles de la carte. Partant de là, il existe des espaces 
disponibles dans ce boîtier dans lesquels on peut disposer la carte sensible incluant 
B,. Nous allons donc fixer pour Xp, Yp et Zp des limites minimales et maximales, 
tenant compte du souci de ne pas « cótoyer » de trop prés la boucle B, et de ne pas 
sortir des limites du boîtier. 

On a ainsi : 

XPmin < [Xp] < XPmax 

Үр < [YP] < Y Pmax 

ZPmin < IZp| < ZP max 

Notre démarche consiste a rechercher pour une combinaison Xp, Yp, Zp (a 
courant donné dans Bp) la valeur de œ autorisant la réduction maximale de la f.e.m. 
induite. 

La figure Ш.36 montre le bornage en Xp, Yp. La variable Zp interviendra plus 
tard en paramétrage de surfaces gauches, lieux des valeurs particulières de œ (notées 
0) pour lesquelles la f.e.m. induite sera minimale. 

A noter que la boucle Bp étant carrée, une symétrie évidente apparait, nous 
permettant de travailler dans un seul quadrant du plan XY. 


Figure 11.36 : Limite de Xp et Yp pour Z < Zpmin 
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Recherche de l'inclinaison optimale 

Sur la base de la représentation de la figure III.36, nous pouvons programmer le 
calcul informatique de manière à obtenir les surfaces gauches, lieux des valeurs de 
Œ, paramétrées par Zp. 

Les figures IIL.37, Ш.38 et IIL.39 donnent (Xp, Yp) dans un quart de plan 
(Xp > 0 et Yp > 0) pour Zp = 2 cm, Zp = 5 cm et Zp = 10 cm. 

Grâce à la mise en œuvre d'un algorithme relativement simple de recherche de о, 
on voit qu'un concepteur peut rapidement trouver l'inclinaison idéale de la carte 
sensible qu'il souhaite inclure dans un boitier de convertisseur. 

Remarque concernant la prise en compte d'un signal i(t) quelconque 

L'étude précédente était basée dans un but de simplification sur un courant de 
fréquence et d'amplitude données. 

Si i(t) est périodique de forme quelconque, il suffit de le décomposer en série de 
Fourier et de raisonner par rang. La courbe 3D de la figure III.40 donne la variation 
de l'amplitude du flux reçu en fonction de œ et de la fréquence. Le calcul a été 
effectué pour le calcul général (bloc 1). On voit que l'amplitude du flux à une valeur 
unique fixée (In) des harmoniques du courant i(t) dans B, ne dépend pas de la 
fréquence, mais exclusivement de a. 

Dans ces conditions l'étude présentée pour la recherche de © est valable pour un 
signal périodique quelconque. 
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Figure 111.37 : Angle optimal œ pour Zp = 2 cm 
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Figure III.38 : Angle optimal о, pour Zp = 5 cm 
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Figures Ш.39 : Angle optimal а, pour Zp = 10 cm 
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Figure Ш.40 : Variation du flux en fonction de la fréquence et de a 


Conclusion du chapitre 


Nous avons vu au chapitre précédent que sous certaines conditions, acceptables 
dans le contexte de notre étude, on peut apporter des simplifications importantes 
dans la démarche de recherche du champ magnétique induit par une boucle de 
câblage d’une carte de convertisseur de puissance. Ces simplifications peuvent aller 
jusqu'à la prise en compte d’une équivalence possible avec une «boucle 
élémentaire » et le choix des équations du champ proche. L'objectif du présent 
chapitre était de montrer comment évaluer les effets de ce champ sur une boucle 
« potentiellement » victime au travers de la recherche de la f.e.m. induite dans celle- 
ci. Dans un souci de rigueur, la recherche de la f.e.m. a été présentée suivant un 
organigramme d'investigation incluant la démarche simplifiée rappelée ci-dessus, 
mais également une démarche générale sans simplifications importantes de manière 
à observer les divergences possibles entre les résultats obtenus. 

On peut conclure que sous réserve de respecter des contraintes de dimensions et 
d'éloignement, autorisant l'étude en champ proche, les équations simplifiées 
permettent d'obtenir la valeur de la f.e.m. avec une erreur raisonnable. 

Aprés avoir calculé ou évalué la f.e.m. dans la boucle « victime », nous avons vu 
comment mener une recherche d'optimisation de la position puis de son orientation 
de maniére à réduire la pénétration du champ et par suite la f.e.m. induite. 

Nous avons ainsi démontré que des zones précises dans l'espace, que nous avons 
qualifiées de « zones tranquilles » pour lesquelles la f.e.m. recueillie sera minimale, 
existent à proximité de la boucle émettrice. 
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Si cela est possible, l'assembleur pourra fixer sa carte sensible au mieux, sinon il 
pourra obtenir par des programmes à complexité réduite, en un point donné de 
l'espace qu'il aura choisi, les conditions d'orientation optimales à prendre en 
compte. 

En résumé nous avons vu comment apporter une aide au positionnement de cartes 
sensibles dans un environnement électromagnétique agressif. 
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CONCLUSION GÉNÉRALE 


Dans cet ouvrage les auteurs ont présenté une contribution à la recherche du 
meilleur positionnement d’un circuit électrique sensible dans l’environnement 
immédiat d’une carte de puissance (convertisseur en 2D). 


L'approche se veut globalement la plus simple possible et répond directement aux 
interrogations issues des tests de mesure en champ proche sur les cartes de 
puissance. 


Par de multiples comparaisons et tests de validité, les auteurs ont prouvé que l'on 
peut adopter une démarche trés simplifiée sous certaines conditions précises à 
respecter. 


L'intérét pour le concepteur est une réponse approximative mais quasi immédiate 
à la question « où placer ma carte sensible ? ». 


Bien entendu le champ magnétique a été essentiellement retenu, car le sujet 
concerne principalement des cartes de puissance, donc des fréquences limitées et des 
niveaux de courant significatifs. 


Ce livre est une contribution à laquelle on pourrait en adjoindre une autre par une 
démarche « duale » concernant alors le champ électrique, et donc les tensions sur les 
pistes ou cábles et les courants induits dans les pistes « victimes » d'un circuit 
sensible. 
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